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USE  OF  SYMBOLIC  MANIPULATION  TECHNIQUES  TO  EXAMINE 
GRAVITATIONAL  AND  UNIFIED  FIELD  THEORIES 

FINAL  TECHNICAL  REPORT 

0.  INTRODUCTION 

This  is  the  final  technical  report  for  work  performed  by  Per¬ 
ception  Technology  Corporation  under  Contract  No.  DAH15  73  C  0369. 
Three  main  objectives  of  the  work  defined  as  "Scope  of  Work"  in  the 
contract,  Item  No.  0001  have  been  satisfactorily  completed.  These 
three  items  are: 

(1)  Verification  of  the  conjecture  that  a  functional  expansion 
(equations  7-9  of  the  proposal)  satisfies,  up  to  a  divergence,  the 
field  equations  (11  of  the  proposal).  Using  MACSYMA  we  have  proven 
that  the  conjecture  is  indeed  correct.  The  proof  is  general  and 
covers  the  full  time-dependent  case  for  the  first  and  second  order 
expansions  of  the  metric.  A  detailed  description  of  this  verifica¬ 
tion  is  given  in  Section  1. 

(2)  Verification  of  the  conjecture  that,  when  conservation  con¬ 
ditions  5^(/-g  ou^uV)  =  0  are  imposed,  the  vanishing  covariant  di¬ 
vergence  of  (11)  reduces  to  the  geodesic  equations  of  motion  (12). 
Using  MACSYMA  we  have  proven  that  this  conjecture  is  also  correct. 


The  computation  covers  the  full  time-dependent  case  and  includes 
first  and  second  order  expansions  as  promised.  Details  of  the 
computation  are  given  in  Section  2. 

(3)  Study  of  a  Unified  Field  Theory  approach  described  in  the 
contract  is  carried  into  first  and  second  order  expansions  of  the 
functional  metric.  The  theory  works  out  satisfactorily  in  the  first 
order  approximation.  In  the  second  order  the  stress-energy  tensor 
is  found  to  differ  from  Maxwell's  tensor  by  an  extra  term.  Vie  have 
not  been  able  to  remove  this  term,  nor  do  we  understand  fully  its 
significance.  Because  of  this  term  the  trace  of  the  stress-energy 
tensor  does  not  vanish.  Hence,  the  Lagrangian  does  rot  reduce  to  its 
special  relativistic  limit  of  correspondence.  However,  we  are  not 
abandoning  the  unified  theory  as  unworkable.  Already  it  comes 
close  to  being  satisfactory  so  that  efforts  toward  a  reasonable  re¬ 
interpretation  is  continuing.  The  details  of  the  computations 
appear  in  Section  3. 

(4)  In  the  course  of  the  above  calculations  the  capabilities 
of  MACSYMA  had  to  be  extended  by  our  compa  v  in  several  respects. 
These  extensions  are  all  important  from  the  point  of  view  of  MACSYMA 
as  a  general  mathematical  tool  and  are  described  in  Section  4. 

(5)  Several  physical  problems  were  investigated  as  additional 
Items  not  listed  in  the  contract.  These  appear  in  Appendices  I-VII 
and  have  to  do  mainly  with  an  overall  comparison  of  the  new  theory 
vis  £  vis  the  conventional  theory  of  Einstein.  The  final  Appendix 
represents  a  concise  formulation  of  the  new  theory  of  gravitation  as 
we  understand  it  at  the  present  time. 
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1.  VERIFICATION  OF  THE  FIELD  EQUATIONS 


The  field  equations  are  given  by 


Rv  -  4 
y  2  p 
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(1.1) 


where  the  left  hand  side  is  the  Einstein  tensor,  the  expression 

□  V  =  ga®4> V |  is  the  general  tensor  d’Alembertian.  The  slash 

y  y  j  op 

denotes  covariant  differentiation  and 
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is  the  stress  energy  tensor  of  the  gravitational  field.  The  quantity 
3^^°  is  a  particular  ordinary  divergence  whose  significance  is  de¬ 
scribed  below.  In  the  theory  the  components  of  g^y  are  functions  of 
the  field  tensor  and  satisfy,  to  second  order,  the  relation 


g  -  n  (1  +  2$  +  2<J>2)  -  4d> 
°pv  pv  r  r  *  ypv 


+  8<pa<l> 
pv  p  w 


(1.3) 


Here  4>  is  the  trace  of  <Ji  and  n  is  the  Lorentz  matrix.  A  general 

v  pv  ° 

relation,  of  which  (1.3)  is  the  second  order  expansion,  is  given  in 
Appendix  VII.  An  additional  equation  in  the  theory  is 


3  <T  =  0 
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This  is  the  "Lorentz  condition"  which  implies  the  local  isotropic 
constancy  of  the  velocity  of  light.  This  physical  condition  is  not 
present  in  the  Einstein  theory  except  for  the  linearized  (first 
order)  Einstein  solution  where  it  appears  as  a  coordinate  condition. 

Our  purpose,  using  MACSYMA  was  to  prove  that  substitution  of 
(1.3)  into  the  left  hand  side  of  (1.1)  results,  to  second  order,  in 
the  right  hand  side  of  (1.1). 

The  definition  of  the  Ricci  tensor  in  terms  of  the  metric  tensor 
is  given  by 


<  *  «"  [»aO  -  ».0  +  <>eS>  - 


(1.5; 


where  the  Christoffel  symbol  {  a}  =  4  g^a(3  g  .  +  3  g  .  -  30g  ). 

J  pa  2  6  '  a6p8  p6og  66pa/ 

The  calculation  proceeds  in  the  following  manner:  The  relation 

(1.3)  is  rewritten  as 


g  =  n  +  X  (24>n  -4*  )  +  X2(2<f>2n  +8<{>a<j>  )  (1.6) 

pv  pv  pv  pv  pv  pv  p  ov 

where  X  is  an  expansion  parameter.  To  compute  Christoffel  symbols 
one  first  needs  the  contravariant  metric  tensor  g^V  which  is  easily 
found  from  (1.6)  via  the  relation 

8pv8W  =  6p  +  °<X3>  <l-7) 


« 
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One  finds 

gWV  =  nPV  -  X(2^WV  -  4*UV)  +  X2(24>2nUV  -  8**PV  +  8^0  (1.8) 

The  Christoffel  symbols  as  well  as  the  Ricci  tensor  may  then  be  com¬ 
puted  in  a  straightforward  manner.  Substitution  of  (1.6)  and  (1.8) 
into  (1.5)  results  in  an  expression  which  contain  terms  of  order  higher 
than  the  second.  Accordingly,  upon  expansion  of  the  Ricci  tensor,  all 

terms  containing  XT1(ri  >  2)  are  dropped.  The  resulting  second  order  ex- 
v 

par.sion  for  is  found  to  have  144  terms.  This  number  would  have 
been  literally  thousands  if  MACSYMA  did  not  have  indicial  manipulation 
capability. 

One  of  the  reasons  for  the  relatively  large  number  of  terms  at 
this  stage  is  the  implementation  of  a  program  to  avoid  the  use  of 
dummy  indices  more  than  twice.  Thus,  each  time  a  dummy  index  occurs 
twice  a  counter  names  the  next  pair  of  indices  with  a  larger  value. 

Once  the  final  expression  is  obtained  a  function  called  rename  cause 
dummy  indices  in  each  element  of  the  expression  to  be  renamed  beginning 
with  the  counter  1.  Thus  many  of  the  terms  cancel  in  the  simplifier 
without  the  implementation  of  any  symmetry  properties.  We  find  rename(R^) 
reduces  the  expression  to  124  terms.  The  next  step  is  the  implementation 
of  the  Lorentz  condition  (1.4).  Many  terms  vanish  as  a  consequence  of 
this  relation  and  we  are  left  with  50.  The  final  step  is  the  contrac¬ 
tion  of  the  terms  with  respect  to  the  Lorentz  matrix.  We  are  then  left 
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with  31  terms  in  the  expression  o  ..  We  form  the  scalar  curvature 

R  -  R^  and  apply  the  functious  Lorentz,  contract  and  rename  once  more. 

We  are  left  with  27  terms  in  R.  The  resulting  expression  for 

RU  -  4  5^R  =  contains  56  terms  as  not  all  terms  are  unique.  This 
v  2  v  v 

expression  of  the  Einstein  tensor  is  further  manipulated  and  reduced 
to 
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(1.9) 


where  □  is  the  special  relativistic  d'Alembertian  defined  by 


□ 


2  = 
0 


(1.10) 


At  this  stage  it  is  convenient  to  check  that  (1.8)  is  correct. 
This  may  be  accomplished  by  computing  the  covariant  divergence  of 
(1.8).  By  definition 


3  GV  -  {  G°  +  {  “}  G° 


y  v 


acr 


va 


(1.11) 
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and  it  is  well  known  that  the  Einstein  tensor,  for  any  metric,  has  a 
vanishing  tensor  divergence.  Using  MACSYMA  we  have  explicitly  proven 
that  divergence  is  indeed  zero. 

If  we  now  examine  (1.9)  in  detail  it  is  seen  that  the  sum  of  the 

first,  second,  third  and  fourth  terms  is  2tV  as  desired.  7t  remains 

P 

to  prove  that  the  remaining  terms  in  (1.9)  are  equal,  within  an  or¬ 
dinary  divergence,  to  the  general  tensor  d'Alembertian  of  <(>^  for  the 
metric  (1.3).  From  the  definition  (1.1)  it  is  found 


□  V  -  +  4<t>“a  aV  -  24>  n^4>v  -  83fi6  aV6 

y  0Yy  Yg  a  y  0  y  B  ay 


26  n  6  +  23  63  6  +  26  C  6  "h  43„6  3  6/> 

yyu0^a  *  u  a  a  (Ty  B  ya  VB 


+  23a<j>  3V<J>  +  43  <t>  D3^aV 

ya  y  aB 


(1.12) 


Substituting  (1.12)  and  (1.2)  into  (1.9)  it  is  found 


g"  -  20 V  -  2t^  =  -  4<j>°363V6  ,  +  83,6,  3%v6  -  43  6 


'ay 


B  ya 


y  aB 


3S*aV 


43  *  5 V8  -  44°363  -  43%3  (iV  -  43%  3%  (1.13) 
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Employing  the  Lorentz  condition  (1.4)  it  is  easily  seen  that  the  right 
hand  side  of  (1.13)  may  be  written  as  a  divergence,  where 
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vo 


. ,Ba.v  q ,  .8,va  .,8V.  .a  .  ~v  Bo 
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"4*B°Vb  -  H*y  -  44»3V*“ 


(1.14) 


Hence,  we  have  established  (1.1)  and  our  conjecture  is  verified.  We 
may  emphasize  again  that  this  computation  would  be  literally  prohi¬ 
bitive  if  carried  out  by  hand  in  its  full  generality. 

In  concluding  this  section  it  will  be  worthwhile  to  point  out 
that  in  the  quasi-static  limit  where  the  heavy  particles  generating 
the  field  are  slow  the  last  term  of  (1.1)  vanishes  on  account  of  the 
Lorentz  condition,  hence,  the  equations  exhibit  an  exact  correspon¬ 
dence  with  a  special  relativistic  tensor  field  theory. 


t 
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2.  EQUATIONS  OF  MOTION  AND  CONSERVATION  LAWS 
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Since  this  is  a  second  order  quantity  it  follows  that 


t  i  =  3  t  *  3  6n„d)  -  23  d>  „  n  <t> 

Pjv  v  p  p*LJ(r  pva6  u0? 


For  any  second  rank  tensor  we  have 


aB 


s;n  *  k  v*  «;>  -  <j>  % 


as  an  identity.  From  Appendix  VII,  equation  (1.1) 


□  2ij)V  =  4ttou  uV 
P  p 


From  (2.6)  and  (2.7)  we  thus  have 


—  fn  2av\ .  =  k 


4ir 
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However,  from  the  definition  of  the  metric  we  have 


o{  0}  u  u^  =  ct3  $ 
PB  a  p 


2auau^  d  <t>  0 

p  o8 


and  from  (2.7) 


□  2<P  =  4tt(J 


Therefore  from  (2.7)  and  (2.10),  (2.8)  becomes 


(2.5) 


(2.6) 


(2.7) 


(2.8) 


(2.9) 


(2.10) 
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‘  %  V/f!  °VB)  -  3/o0*  +  2V«A*“*  <2-ii) 


Thus,  from  (2.5)  and  (2.11)  we  find  (2.2)  becomes 


3  (/g  au  uV)  +3  3  eaB  =  o 
v  W  a  6  y 


(2.12) 


oince  is  also  a  second  order  quantity.  This  is  the  conservation 

law  for  the  new  theory.  In  the  quasi-static  limit  3  3  e°®  =  Gftf3-) 

a  6  u  ' 

and  «e  ara  left  vich  3v(/g  duuV)  -  0.  He  believe  that  for  non-static 
metrics  3^^  is  involved  in  gravitational  radiation.  This  subject 
will  be  investigated  in  the  next  phase  of  the  contract. 


Let  us  now  return  to  (2.2)  which,  from  (2.8)  yields  (to  second 


order) 


8"  (k  V/g  OU/>  -  °  U>  “8“a)  +  +  VbS*  (2.13) 


From  ('’.12)  this  reduces  to 


°  U}  u^ua  +  96  "  0 


(2.14) 


However,  from  correspondence  with  special  relativity  and  other  local 
field  theories  the  equations  of  motion  of  test  particles  are  given  by 


a  ^7  =  3v  <*>> 


(2.15) 
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Hence,  from  (2.14)  and  (2.15)  we  find  the  geodesic  equations  of 
motion 


d2x 
_ y 

ds2 


r  a,  e 

V  "  “a 


(2.16) 


which  are,  of  course,  equivalent  to  their  more  familiar  form 


d2xy 


ds 


.  t  y-.  a  6 
+  t  u  u 
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(2.17) 


Thus,  we  have  verified  that  the  geodesic  equation-  of  motion  (2.17) 
follow,  up  to  a  divergence,  from  the  field  equations  (2.2)  upon  ap¬ 
plication  of  the  conservation  laws  3v(/-g  ou^uV)  =  0.  We  have  thus 
verified  the  conjecture  that  in  first  and  second  order  the  geodesic 
equations  of  motion  are  a  consequence  of  the  field  equations  plus 
the  conservation  laws.  This  statement  is  rigorously  valid  in  the 
quasistatic  limit.  In  the  most  general  fully  time-dependent  case 
the  divergence  terms  seem  to  be  related  to  the  radiation  reaction 
and  radiative  damping.  In  the  quasi-static  limit  radiation  tends 
to  zero,  divergence  terms  vanish,  hence  the  geodesic  equations  of 
motion  become  rigorously  valid. 
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3.  SKEW-SYMMETRIC  $  AND  UNIFIED 

yv 

FIELD  THEORY 


Our  methods  in  the  third  phase  of  the  contract  do  not  follow 
the  standard  approach  to  unified  field  theories.  Our  philosophy 
is  as  follows:  We  know  that  to  second  order  (Appendix  7,  (2.6)), 
the  definition 


g 


n 


e2($  -  2*) 


(3.1) 


produces  a  viable  theory  of  gravitation  when  the  Einstein  tensor  G 
is  evaluated  for  symmetric  The  obvious  question  to  be  considered 

concerns  the  form  G  takes  when  is  skew-symmetric.  Skew-symmetry 
of  <f>^  implies  <J>  =  0,  hence  we  consider  the  metric  tensor  to  take  the 
form  (to  second  order) 
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yv 


n  +  Ae 
yv  yv 


+  4A2 


(3.2) 


where  n  ,  is  the  symmetric  Lorentz  matrix  but  now  e  =  -  e  The 
yv  yv  vy 

parameter  a  in  (3.2)  takes  the  value  +  1  and  can  only  be  determined 
later  via  a  correspondence  principle.  A  is  an  expansion  parameter 
such  that  A  <<  1.  We  do  not  require  the  second  order  contravariant 
metric  tensor  but  will  need  the  first  order  part.  It  is  found  that 


g 


yv 


„yv 


-  Ae 


yv 


(3.3) 
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satisfies  guvgVp  =  6P  +  0(1*)  and  hence  (3.3)  is  the  relation  we 
adopt. 

The  only  postulate  we  shall  employ  is  similar  to  other  unified 
theories.  Namely,  the  metric  tensor  is  to  have  a  vanishing  covariant 
derivative.  Thus 

8uv;a  "  8a8yv  “  8avrua  “  8yarav  =  0  (3.4) 

From  this  relation  we  may  now  obtain  the  connection  coefficients  ra 

ua 

for  the  metric  (3.2)  which  we  write  in  the  more  convenient  form 


!„v  -  V  +  X£„v  +  X\ 


uv 


(3.5) 


Since,  ir.  metrical  theories,  connection  coefficients  are  linear  in 
first  derivatives  of  the  metric  and  since  is  constant  with  respect 
to  coordinate  differentiation  it  follows  that  to  first  order  (3.4) 
and  (3.5)  imply 


3  e 
a  yv 


-  n  r"  -  n  r  =  o. 

av  ya  ya  av 


The  solution  of  this  relation  is 


(PP»<0 


(3.6) 


(3.7) 
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where  we  adopt  the  notation 


(yp.a)e 


1 

2 


(3  e  +  3  e 
a  pp  y  pa 


3  e  ) 
P  o\i 


(3.8) 


Using  an  identical  method,  with  the  notation  in  (3.8)  we  find  the 
second  order  expression  for  to  be 


rpa  =  nP°  l^p>a)e  + 


(3.9) 


where 


3  ip 

a  UP 


=  3  h 
a  up 


(u0,a)E 


(ae,p)e  . 


(3.10) 


Substitution  of  (3.10)  into  (3.9)  results  in  25  terms  for  ra  .  This 

ya 

expression  is  all  we  need  to  examine  geometrical  tensors  in  this 
framework.  This  is  as  far  as  our  calculation  proceeded  by  hand. 
MACSYMA  was  then  given  the  definition  (3.9)  and  (3.10). 

We  first  note  that  the  first  set  of  Maxwell's  equations 
can  be  written  as 


3v(»/-g  gyV)  -  4,tJV  (3.11) 

Then  the  second  set  is  obtainable  as 


3o8By  +  3S8Y*  +  \8cB 


0 


(3.12) 
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From  (3.8)  one  can  further  prove  the  Lorentz  force 


d2x 
_ P 

ds2 


{  a\  e 
l  0i  u  u 
pB  a 


F  J 

pa 


(3.13) 


One  can  therefore  see  that  the  first  order  equations  of  electrody¬ 
namics  can  be  accommodated.  The  problem  then  reduces  to  a  second 
order  calculation  to  see  if  the  Maxwell  tensor  is  also  obtainable. 

Given  r°  to  second  order  we  shall  construct  the  Ricci  tensor 
pa 

using  MACSYMA,  via  the  definition 


rs  ,  ?s  „t 
is,k  ts  ik 


rs  r1" 

tk  is 


(3.14) 


In  the  standard  version  of  Tensor  Fasl  (Section  4),  (3.14)  may  be  ex¬ 
panded  to  second  order  with  no  difficulty  as  MACSYMA  does  not  sort 
indices.  The  expression  for  R^  contains  102  terms  which  is  reduced 
to  90  using  rename.  A  new  version  of  Tensor  Fasl  was  written  to 
enable  the  contraction  of  expressions  when  the  order  of  indices  is 

important  as  in  this  case.  Contracting  R^,  employing  the  Lorentz 
y  v 

condition  e  ^  =  0  and  grouping  all  divergence  we  find  that 


R  "7  [3  e  . 
pv  4  p  aB 


3  eaB  -23  e  36e  a  -  2Q*e  ]  +  Div.  (3.15) 
v  6  pa  v*  0  pv 


2  _  UV 

where  □  q  =  n  3^3^.  This  result  has  been  obtained  by  requiring  that 
the  expression  for  R^y  reduce  to  the  symmetric  case  ($  =  0)  when 
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e  *  "  Evu*  ^*is  t^'e  correspondence  limit  which  gives  a  *  1 

in  (3.2). 

From  (3.12)  and  (3.3)  it  follows  that 


<nyv  -  xeyv)R 


yv 


-4x2£ 


_yv,o 


yv,a 


(3.16) 


From  (3.12)  and  (3.13)  we  find 


G  =  y  [3  e  3  e013  -  23fte  3Be  a  -  2D*e  -  \  n  3  e  Q3ae“B]  +  Div. 
pv  4  p  ag  v  p  pa  v*  0  pv  2  pv  a  a0 


(3.17) 


Note  that  this  is  precisely  the  same  result  we  obtain  in  the  symmetric 

case  if  we  set  4>  =  0  and  choose  e  *  -  44 

yv  vp 

Note  that  the  replacement  of  D2  instead  of  d2  would  only  bring 

t » 

antisymmetric  terms  hence  would  not  alter  the  essential  character  of 
the  result.  Two  of  the  three  non-linear  terms  in  (3.14)  are  exactly 
of  the  form  of  a  Maxwell  stress-tensor.  The  only  term  which  is 
troublesome  is  the  first  one.  At  the  present  time  we  do  not  know  what 
to  do  with  this  result  except  to  say  that  the  whole  scheme  almost  works 
out  but  doesn't  quite  make  it  because  of  the  extra  term.  Future  in¬ 
vestigations  may  indicate  a  way  to  overccme  this  problem  or  help  us 
to  clarify  the  situation. 


In  conclusion  although  the  unified  field  theory  direction  did 
not  yet  work  out  to  our  satisfaction  it  showed  sufficient  promise 
that  it  may  do  so  by  a  suitable  reinterpretation.  In  any  case  it 
will  be  fair  to  say  that  with  the  help  of  the  M..CSYMA  we  were  able 

to  carry  the  subject  further  than  any  other  calculation  made  in  the 
past. 


4.  NEW  CAPABILITIES  OF  MACSYMA 
THE  TENSOR  PACKAGE 


The  MACSYMA  system  consists  of  over  250  functions  for  the  mani¬ 
pulation  of  mathematical  expressions.  There  are  not  only  functions 
for  performing  mathematical  computations  but  also  those  for  many 
necessary  miscellaneous  manipulations  such  as  extracting  parts,  sim¬ 
plifying,  input  output,  etc.  Each  one  of  these  if  made  up  of  one  or 
as  many  as  several  hundred  programs  (in  the  case  of  integration) 
written  in  the  LISP  language.  The  expressions  dealt  with  are  composed 
of  operators  and  arguments  to  them  which  are  numbers,  variables,  func¬ 
tions  references,  arrays,  lists,  and  matrices.  The  mathematical 
knowledge  built  for  MACSYMA  is  continually  expanding  as  more  and 
more  programs  are  written  to  provide  additional  capabilities.  At 
present  there  are  over  3000. 

If  one  wishes  to  introduce  a  new  mathematical  entity  into  the 
system  then  he  would  not  only  have  to  cast  it  in  terms  of  the  existing 
syntax  but  also  it  would  be  necessary  to  provide  the  programs  to  per¬ 
form  the  desired  manipulations  on  it  as  it  would  be  foreign  to  most 
of  the  existing  programs.  Perception  Technology  Corporation  has 
embarked  on  an  effort  to  create  facility  for  dealing  with  tensor 
manipulation  as  MACSYMA  had  nothing  that  would  provide  a  convenient. 


efficient,  and  workable  mechanism  for  indicial  types  of  manipulations. 
The  basic  program  we  have  created  is  called  "Tensor  Fasl".  This 
program  is  now  quite  sophisticated  and  we  would  estimate  it  contains 
40%  of  the  tensor  analysis  found  in  standard  texts.  However,  it  is 
capable  of  solving  problems  in  tensor  calculus  that  are  impractical 
to  do  by  hand. 

Since  tensors  are  defined  not  by  the  way  they  are  written  but  by 
certain  properties  they  have,  we  represent  them  as  more  general  quanti¬ 
ties  which  we  call  indexed  objects.  They  possess  a  name  and  two  lists 
of  indices,  the  first  list  being  the  subscripts  (covariant  indices) 
and  che  second  list  being  the  superscripts  (contravariant  indices) , 
These  indexed  objects  are  represented  as  functions  of  two  arguments, 

which  are  the  above  mentioned  lists.  If  either  set  of  indices  is 

j  ilk 

absent  it  is  represented  by  an  empty  list.  Thus  6^  and  A  would  be 
represented  as  5([i],[j])  and  A([  ],[i,j,k]). 

Ordinary  differentiation  of  an  indexed  object  (say  E)  with 

k 

respect  to  the  coordinate  x  is  obtained  by  the  command  DIFF(E,k), 

Since  MACSYMA  wouldn't  know  that  E  depends  on  x  it  would  give  0  but 
the  differentiation  program  in  Tensor  Fasl  has  been  modified  to 
assume  that  all  indexed  objects  depend  on  all  coordinates.  If  an 
indexed  object  is  independent  of  all  coordinates  this  is  stated  by 
the  command  DECLARE (E, constant)  and  subsequent  differentiation  of  E 
would  yield  0.  The  derivative  of  an  indexed  object  causes  the 
coordinate  with  respect  to  which  differentiation  is  carried  out  to 


be  appended  to  the  function  which  represents  the  indexed  object  in 
MACSYMA.  The  indices  which  denote  differentiation  with  respect  to 
a  coordinate  are  sorted  in  alphabetic  order  to  take  advantage  of 
the  commutativity  of  differentiation  in  order  to  simplify  expressions 
composed  of  those  indexed  objects. 

At  present  three  basic  features  have  been  implemented  -  symmetries 
and  contractions,  and  covariant  differentiation.  For  the  first  three 
functions  are  provided.  One  for  declaring  symmetries  (e.g.  G({  ],[j,i]) 
is  equal  to  G([  ],[ijj]))>  one  for  removing  them,  and  one  for  dis¬ 
playing  them  in  case  the  user  forgot  which  ones  he  declared.  For  the 
purpose  of  contracting  products  of  indexed  objects  in  expressions, 
four  functions  were  provided.  One  to  declare  that  some  indexed  object 
contracts  with  another  to  form  a  third,  another  function  to  remove 
these  declarations,  another  to  display  them,  and  also  a  function  to 
perform  the  contractions.  This  latter  function  is  rather  difficult 
to  code  as  it  is  necessary  to  make  several  passes  over  the  expression 
to  make  sure  that  all  possible  contractions  have  been  utilized.  The 
function  also  performs  the  substitutions  implied  by  the  Kronecker 
delta  function  if  it  occurs  in  an  expression.  In  addition  a  function 
'Show'  has  been  written  to  display  expressions  containing  indexed 
objects  in  more  natural  notation  than  MACSYMA  ordinarily  would. 

MACSYMA  assumes  that  A([i,j],[k],l,m)  is  a  function  of  four  arguments, 
and  would  display  it  as  such,  whereas  we  like  to  see  it  displayed 


as  Tensor  Fasl  also  performs  multiple  covariant  differentia¬ 

tion  on  indexed  objects  with  any  number  of  indices.  These  expressions 
may  be  expanded  in  terms  of  Christoffel  symbols  of  the  second  kind 
(and  their  derivatives).  In  addition,  numerical  tensors  are  defined 
as  well  as  the  relation  between  connection  and  metric. 

The  following  list  of  capabilities  of  the  MACSYMA  system  is  di¬ 
rectly  attributable  to  the  combined  effort  of  Perception  Technology 
Corporation  and  our  program  consultants  at  Project  MAC.  At  the 
start  of  our  contract  the  MACSYMA  system  was  only  capable  of  computing, 
for  a  given  metric,  the  inverse  metric,  Christoffel  symbols  and  Ricci 
tensors.  The  latter  of  these  had  a  bug  which  was  discovered  and  cor¬ 
rected  by  Perception  Technology  Corporation.  The  list  of  capabilities 
as  of  June  15,  1974  include: 

1)  For  any  metric  tensor  the  Ricci  tensers  and  Einstein  tensors 
(covariant  or  mixed)  can  be  computed  in  CRE  (canonically  rational  ex¬ 
pression)  form. 

2)  For  any  metric  tensor,  the  inverse  metric,  Christoffel  symbols 
of  both  kindr,  and  mixed  Einstein  tensors  can  be  computed  in  CRE  or 

in  a  truncated  power  series  to  arbitrary  order  in  an  expansion  param¬ 
eter.  We  are  limited  only  in  the  memory  of  the  system. 

3)  For  any  metric  tensor  the  components  of  the  full  Riemann 
tensor  can  be  computed  in  CRE  or  in  a  truncated  power  series. 
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4)  For  any  metric  tensor  the  components  of  the  Weyl  tensor 
can  be  computed  in  CRE  or  in  a  truncated  power  series.  We  have 
not  attempted  classifications  yet  although  this  facility  is  now 
not  difficult. 

5)  For  any  metric  tensor  the  equations  of  motion  can  be  com¬ 
puted  in  CRE  or  in  a  truncated  power  series. 

6)  For  any  metric  tensor  the  d ’Alembertian  and  the  of  the 

y 

new  theory,  both  applied  to  a  scalar  field,  can  be  computed. 

7)  For  any  metric  tensor  the  components  of  any  second  rank 
tensor  in  a  transformed  coordinate  system  can  be  computed. 

8)  For  a  function  $(z-t)  one  has  the  differential  identities 

$  -  “  <t>  »  <j>  „  =  -  <J>  ,  <i>  =  +  A  _  which  are  the  wave  conditions 

z  t  zt  tt  zz  tt 

of  certain  specially  chosen  metrics.  For  components  of  the  metric 
tensor  containing  functions  of  this  kind  the  wave  conditions  may 
now  be  automatically  implemented  with  MACSYMA  thus  tremendously 
shortening  computations  concerned  with  gravitational  radiation. 

9)  Given  a  metric  tensor  in  the  form  g  =  n  +  +  X2i|i 

°yv  yv  yv  Tyv 

where  is  the  Lorentz  matrix,  X  is  an  expansion  parameter,  and 

$  and  \p  are  arbitrary  tensor  fields,  MACSYMA  is  now  able  to  com- 
yv  yv  J 

pute  the  Riemann,  Ricci  and  Einstein  tensors  and  their  invariants 
to  second  order  in  X.  This  last  item  was  especially  important  for 
the  completion  of  the  general  calculation  and  represents  the  only 
existing  tensor  manipulation  package  according  to  computer  experts. 


5.  APPENDICES 


Below  are  included  seven  appendices  dealing  with  subjects  in¬ 
vestigated  during  the  contract  period.  While  these  are  not  strictly 
of  the  identical  subjects  as  outlined  in  the  contract,  they  represent 
advances  in  the  understanding  of  the  new  theory  as  well  as  Einsteinian 
theory.  Appendices  I,  II  and  V  deal  with  consequences  of  the  general 
spherically  static  solution  of  Einstein's  vacuum  equations,  and  each 
has  been  submitted  for  publication.  Appendix  VI  gives  the  general 
static  solution  of  the  vacuum  equations  for  a  metric  with  plane  sym¬ 
metry.  Appendix  III  deals  with  Plane  Gravitational  Waves  and  repre¬ 
sents  an  extension  of  a  paper  read  at  the  American  Physical  Society 
Meeting,  April  1974.  Appendix  IV  displays  some  of  the  interesting 
properties  of  exponential  matrices.  Finally,  Appendix  VII  represents 
the  culmination  of  our  current  understanding  of  the  new  theory. 

Efforts  will  be  made  to  publish  as  many  as  possible  of  these  results 
and  appendices  in  recognized  scientific  journals  and  conference 
proceedings. 


APPENDIX  I 


GENERAL  SOLUTION  OF  EINSTEIN'S 
EQUATIONS  WITH  SPHERICAL  SYMMETRY* 


ffllseyin  Yilmaz 

Perception  Technology  Corporation,  Winchester ,  Massachusetts  01890 


ABSTRACT 

A  general  solution  for  Einstein's  theory  of  gravitation  corres¬ 
ponding  to  spherical  symmetry  is  presented .  It  is  shown  that  the 
usual  Schwarzchild  metrics  are  members  of  this  general  solution. 

Under  the  known  computational  methods  believed  to  be  generally  cova¬ 
riant  some  of  the  member  solutions  to  be  exhibited  seem  to  make  little 
or  no  physical  sense.  One  that  seems  to  be  mathematically  quite 
unique  and  interesting  for  possible  astrophysical  applications  is 
described  in  detail. 

By  the  use  of  a  well-known  method  of  tensor  transformation  it  is 
proven  that  the  line  element 

ds 2  =  fc2dt2  -  (jZfJ  (dQ2  +  sin2Qd$ \2)  -  f  “p )  dr2  (1) 

where  f  is  a  thrice  differentiable  function  of  r  is  a  general  solution 
of  Einstein' s  field  equations  R^  -  0.  The  usual  non-isotropic  and 
isotropic  Schwarzchild  solutions 1  are  members  of  this  general  solution 
with  the  choices 

*This  research  is  supported  by  the  Advanced  Research  Projects  Agency 
of  the  Department  of  Defense  Under  Contract  No.  DAHC15  73  C  0369. 


f  -  I  -  2* 


f  -  i  -  2*  fi  +  $/2; 


(2) 


-2 
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where  <p  =  m/r  so  that  for  large  r  the  first  order  Newtonian  corres¬ 
pondence  is  satisfied.  Similarly  f  =  (1  -  $)/ (1  +  <f>)  yields  the 
well-known  Fock  solution 2  satisfying  the  correspondence  condition  in 
first  order.  The  solution  (1)  is  new  and  offers  enormous  flexibility 
as  we  shall  see  below.  (It  was  originally  conjectured  during  the  per¬ 
formance  of  an  ARPA  contract  directed  to  study  the  solutions  of  au¬ 
thor'  s  new  theory  of  gravitation^ ,  and  later  explicitly  verified  by 
the  use  of  the  MACSYMA  system  of  Mathematic's  Laboratory ,  Project  MAC, 
M.I.T.) 

Actually  it  is  somewhat  disturbing  to  find  that  f  can  be  so 
general  as  in  (1) ,  because  we  seem  to  be  able  to  set 

f  =  Z  a  < t>n  (3) 

n 
n 

and  produce  an  arbitrary  number  of  Schwarzchild  type  singularities , 
or  choose. 


f  = 


(4) 


and  obtain  a  solution  with  no  singularity  at  all,  except  at  r  =  0. 
These  facts  seem  to  put  a  question  mark  on  the  problem  of  black 
holes.  Previously  we  were  under  the  impression  that  only  one  legiti¬ 
mate  Schwarzchild  surface  can  occur.  The  solutions  corresponding  to 
(3)  contradict  this  contention. 


Note  that  some  of  the  possibilities  contained  in  (3)  might  also 
be  problematic  for  either  the  experimental  or  the  computational  as¬ 
pects  of  the  theory.  For  example,  by  setting 


f  -  1  -  2$  +  Jb$2 


(5) 


ds2  -  (1  -  2$  +  b$2)dt2  -  (1  -  bW2)~2  (r2dQ2  +  r2sin2Qd$2) 


L  n  -  M)2  dr2 

f  (1  -  bW2)h 


(6) 


where  b  is  arbitrary ,  one  can  prove  by  the  use  of  a  standard  computa¬ 
tional  method 4  (which  is  believed  to  be  universally  applicable)  that 
the  three  crucial  tests  of  red  shift,  light  bending  and  perihelion 
advance  are  satisfied,  but  the  radar-echo  delay  depends  on  b  as 


0  0  I 

ht  =  4M  (In  4  —  +  b/2)  (7) 

H2 

Since  b  is  arbitrary  one  may,  for  instance,  set  b  =  -22  and  get  a  b- 
dependent  part  which  cancels  the  whole  gravitational  effect  expected 
from  the  theory.  Of  course  the  experiment  is  more  sensitive  to  the 
differential  delay,  3 t/%R  =  -8M/R,  but  b  can  always  be  set  arbitrarily 
large  to  counter  arguments  based  on  the  relatively  small  sensitivity 
to  b.  (Units  are  c  =  G  =  1. )  This  ambiguity  is,  in  fact,  a  severer 
form  of  a  well  known  difficulty  already  present  between  the  two 


Schwarzchild  solutions  b  =  0  (standard)  and  b  =  2  (isotropic)  given 


by  (2).  As  early  as  1959  the  present  author  insisted  that  the  two 
Schwarzchild  solutions  might  not  be  physically  equivalent,  and 
proposed  a  test  of  the  isotropy  of  space  using  masers.5  In  1966 
Ross  and  Schiff  exhibited  the  difficulty  more  explicitly  and  offered 
a  possible  resolution  by  bringing  in  the  dynamics  of  the  solar  system 
to  bear  on  the  problem .6  But  the  paradox  was  really  never  resolved 
because  the  Ross-Schiff  argument  depended  on  choosing  the  standard 
Schwarzchild  metric  as  the  starting  point  of  calculation.  If  one 
starts  with  the  Non-isotropic  Schwarschild  line  element  and  apply  the 
same  procedure  one  can  show  that  the  original  discrepancy  remains 3, 
namely,  the  isotropic  line  element  still  yields  a  time-delay  20  ys 
larger  than  the  standard  line  element.  In  (6)  this  difference  is 
now  made  to  depend  on  an  arbitrary  parameter  b,  such  that  the  dis¬ 
crepancy  is  St  =  2M6b  hence  it  can  be  arbitrarily  large.  A  recent 
investigation7  has  further  revealed  that  the  choice  f  =  1  -  2$  +  2a<p2ln<p 
leads,  upon  the  application  of  the  general  Robertson-Noonan  computa¬ 
tional  methods1*,  to  correct  first  order  resultr  for  the  red  shift, 
light  bending  and  the  3t/3i?  =  -  8M/R  part  of  the  time-delay  but  gives 
a -dependent  differences  for  the  perihelion  motion  and  the  absolute 
value  of  the  gravitational  time-delay  by  the  factors  1  +  ol/3  and 
1  -  3a/ 2  respectively.  In  view  of  the  claimed  covariance  of  the  com¬ 
putational  method  used  this  result  is  very  surprising  indeed. 

This  strange  situation  may  be  summed  up  as  follows:  "If  general 
relativity  is  correct  all  metrics  (1)  must  be  observationally  equivalent 


and  any  observation  made  using  different  frames  of  reference  (1) 
must  be  capable  of  correlating  without  discrepancy ."  The  above  dis¬ 
crepancies  may  then  be  interpreted  in  two  ways:  a)  We  possess  a 
generally  covariant  calculational  method  but  general  relativity  is 
not  correct  beyond  first  order,  b)  The  theory  is  correct  but  we  do 
not  possess  a  generally  covariant  method  of  calculation.  The  present 
writer  tends  to  side  with  the  first  interpretation.  In  two  recent 
articles 8  it  is  shown  that  general  relativity  conserves  rest-mass 
hence  it  is  generally  covariant  under  a  group  of  transformation  (1) 
that  conserves  rest-mass.  This,  however,  conflicts  with  the  cor¬ 
respondence  limit  of  special  relativity  which  conserves  energy- 
momentum  and  not  rest-mass.  Methods  of  computation  do,  however, 
implicitly  imitate  special  relativity  hence  presume  conservation 
of  energy -momentum.  It  appears  therefore  that  the  clash  between  the 
two  is  reflected  in  the  computed  discrepencies  above. 

From  among  the  new  solutions  contained  in  (1)  the  possibility 
of  setting  or  requiring  f  =  e  ^  is  probably  the  most  interesting 
for  general  relativity .  In  this  case  we  have  the  line  element 

2  -2< |> 

ds2  .  e~2idt2  -  I  — file2  +  sin2ddi2) - — '  8,,-  r  dr2  IS) 

1-e'2*  r*ri-e'2V 

This  new  line  element  has  the  following  remarkable  properties: 
a)  For  r  »  m  it  reduces  to  the  usual  isotropic  Schwarzchild  line 
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element  for  all  experimental  purposes,  hence  it  reproduces  all  four 
of  the  experimental  tests  [in  first  and  second  order  (8)  is  the  same 
as  the  non-isotropic  Schwarzchild  line  element ]  b)  In  the  limit  of 
slow  motion  of  the  test  body  the  geodesic  equations  of  motion  reduce  to 


d2x 

_ y 

dt2 


.  o,  o 
=  {.  }  u  u 
yo  o 


(9) 


which  are  exactly  the  Newtonian  equations  in  spherical  coordinates, 
c)  In  the  limit  of  slow  motion  of  the  source  the  field  equations  reduce 
to 


V2<j>  =  -4  no  (10) 

which  is  exactly  the  Poisson  equation.  By  contrast  ether  line  elements 
reduce  to  these  Newtonian  limits  in  an  approximate  sense  only,  d)  It 
is  the  only  solution  of  Einstein's  equations  which  apply  both  to  the 
exterior  and  the  interior  regions  of  a  mass  distribution  without  a 
change  in  form,  e)  It  is  harmonic  in  the  sense  that  if  the  metric  is 
written  in  the  form,  ds *•  =  Dc2dt2  -  B(r2dQ2  +  r2sin2Bd$2)  -  Adr2 ,  that 
is  as  a  factor  to  special  relativity® ,  then  %r('/-g  g^  )  -  ^r(^DB2A  A  *) 
[in  Cartesian  coordinates  x,  y,  z  (isotropic)  the  line  element  is  not 
harmonic ].  f)  It  has  no  singularity  anywhere  except  at  r  0  where  it 
reduces  to 


ds 2  =  ~(2m)2  (dQ2  +  sin26d<t>2) 


(ID 
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Tiie  latter  represents  the  surface  of  an  inpenetrable  sphere  of  radius 
2m.  This  spherical  surface  seems  to  be  suspiciously  similar  to  the 
Schwarzchild  surface  but  note  the  difference  that  (8)  has  an  essen¬ 
tial  singularity  at  the  surface ,  not  an  ordinary  one.  Furthermore , 
there  is  really  no  interior  of  this  sphere  as  r  starts  out  at  the 
surface.  A  possible  objection  to  (8)  might  be  that,  like  all  other 

solutions  of  Einstein's  theory,  it  excludes  the  vacuum  gravitational 

v 

stress-energy  tensor  t  and  thereby  violates  the  Newtonian  corres¬ 
pondence  because  we  shall  be  missing  the  stress-energy  tensor 


tV  =  - 9  <f>  3V<J>  +  i  6V  3X4>  M  (12) 

p  yr  r  2  u  X 

of  the  Newtonian  theory. 

It  is  a  pleasure  to  thank  Dr.  Richard  Pavelle  and  Mr.  Alphonsus 
Fenelly  for  valuable  discussions. 


*This  research  was  supported  by  the  Advanced  Research  Projects 
Agency  of  the  Department  of  Defense  under  Contract  No.  DAHC15  73  C  0369. 
That  the  conjectural  line  element  (1)  is  a  general  curved-space  solu¬ 
tion  of  Einstein's  equations  was  proven  by  a  tensor  transformation  on 
the  metric.  There  are  also  purely  flat-space  solutions  of  the  form 
ds 2  =  c2dt2  -  F26302  +  sin20d<j>2J  -  F,2dr2  which  are  not  discussed  in 
this  paper.  These  solutions  are  checked  by  MACSYMA,  a  symbolic  mani¬ 
pulation  system,  at  Project  MAC,  MIT.  The  author  is  grateful  to  Dr. 
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cation  possible.  A  direct  proof  of  (1)  is  found  by  Dr.  Richard  Pavelle 
of  Perception  Technology  Corporation  by  using  the  extensive  manipulative 
power  of  MACSYMA  system.  Whether  these  are  the  only  possible  general 
solutions  of  the  theory  in  spherical  coordinates  has  not  yet  been  es¬ 
tablished,  although  a  reasoning  based  on  Birkhoff's  theorem  seems  to 
indicate  that  there  are  no  ethers. 
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APPENDIX  II 

DOES  GENERAL  RELATIVITY  PREDICT  BLACK  HOLES?* 


Richard  Pavelle  and  Huscyin  Yilmaz 
Perception  Technology  Corporation ,  Winchester ,  Massachusetts 


ABSTRACT 


The  general  static  spherically  symmetric  so¬ 
lution  of  Einstein's  vacuum  equations  is  given.  It 
is  shown  that  besides  the  non-isotropic  Schwarzchild 
solution  which  possesses  a  singularity  usually  as¬ 
sumed  to  represent  a  black  hole,  we  may  construct 
equally  valid  solutions  with  a )  any  number  of  black 
holes,  b)  a  black  hole  of  arbitrary  radius,  c)  no 
black  holes.  All  solutions  given  in  this  article 
except  one  satisfy  the  three  experimental  tests. 

The  arguments  presented  herein  question  whether 
General  Relativity  does,  in  fact,  unambiguously 
predict  the  existence  of  black  holes. 


01890 


*This  research  is  supported  by  the  Advanced  Research  Projects  Agency 
of  the  Department  of  Defense  Under  Contract  No.  DAHC15  73  C  0369. 


It  is  known  that  the  general  static  and  spherically  symmetric 
line  element  in  general  relativity  is  of  the  form 1 

ds 2  =  e°  dt2  -  e^r2(d62  +  sin28d$2)  -  eXdr2  (1) 

where  V,  p,  X  are  functions  of  r.  The  vacuum  field  equations ,  =  0, 

of  Einstein's  theory,  however,  yields  three  equations  only  two  of  which 

9 

are  independent. 

*11  -  +  i  v"  +  f  **'  -  7  x'  +  j  u'2  -  J  x'v'  -  7  x'v' 

+  j  V'2  =  0  (2) 

*44-eV“X  (±v"  +jvV  -  ^v'2j  »  0  (3) 

It  is  therefore  clear  that  (1)  must  contain  an  arbitrary  function,  f, 

of  r.  Since  a  general  solution  containing  such  an  arbitrary  function 

is  never  written  down  one  usually  imposes  coordinate  conditions  such 

i  uv 

as  p  =  0  (standard) ,  p  -  X  =  0  (isotropic) ,  or  d^(v-g  g  )  =  0  (har¬ 
monic)  and  obtains  special  line  elements .  All  the  solutions  so  far 
obtained  by  such  specializations  have  a  common  property  in  that,  by 
a  method  of  calculation  believed  to  be  applicable  to  any  line  element 
with  arbitrary  interpretations  of  t  and  r,  they  yield  the  three  crucial 
tests  of  red  shift,  light  bending  and  perihelion  advance  within  experi¬ 
mental  accuracy2 . 
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We  have  recently  determined 3  the  general  solution  of  the  diffe¬ 
rential  equations  (2)  and  (3).  The  form  of  the  line  element  (1)  con¬ 
taining  an  arbitrary  function ,  f ,  is  found  to  be 

ds2  =  fdt2  -  (  )2(d&2  +  sin2Qd$2)  -  j  [  (  )'  dr2  (4) 

This  general  solution,  v/hich  can  be  checked  by  using  equations  (2)  and 

(3)  yields  the  above  specializations  directly  by  specifying  f  as 

-2  -1 

f  =  1  -  2 <J>,  f  =  1  -  2§(1  +  §/2)  ,  f  =  1  -  $(1  +  <j>J  where  <J>  =  m/r. 

What  is  found,  by  a  general  application  of  (4)  is,  however,  quite  dis¬ 
turbing.  We  shall  first  give  a  few  examples  to  indicate  the  nature  of 
our  worries  and  later  comment  on  the  possible  implications  of  the  si¬ 
tuation  for  the  theory  as  a  whole. 

i)  All  the  special  solutions  so  far  knovm  involve  only  one  singu¬ 
larity  of  the  Schwarzchild  type.  But  consider  the  choices 

P  P 

eV  =  f  =  n  (1  -  C  4)  ,  l  C.  =  2  ( 5 ) 

3  -  1  J  3  3 

where  C >  0.  These  solutions  are  exact,  satisfy  the  boundary  conditions 
at  infinity,  reduce  to  the  Newtonian  theory  in  first  order  and  yield  the 
three  crucial  tests  correctly  but  in  general  possess  p  (an  arbitrary  number) 
singularities.  If  one  computes  the  location  of  the  null  hypersurfaces 
corresponding  to  these  singularities 3  it  is  found  that  these  singularities 
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satisfy  the  usual  definition  of  black  holes.  However,  we  can  have 
p  such  singularities,  p  being  any  number.  Which  one  of  these  singular 
surfaces  will  be  identified  with  the  physical  black  hole  assumed  to 
exist  in  the  literature? 

ii)  Another  choice  which  satisfies  che  above  desiderata 

eV  =  f  =  (1  -  2W1/K  (6) 

where  £  is  a  positive  number.  In  this  case  we  have  only  one  singularity 
which  corresponds  to  a  null  hypersurface  but  its  position 

r  =  2mE  (7) 

can  be  anywhere.  For  E  -  2.5  I05  the  singularity  will  engulf  the 
sun  without  affecting  the  three  crucial  tests  theoretically .  Is  the 
sun  a  black  hole? 

Hi)  Consider  now  E  to  be  negative,  namely, 

v  -7/E 

e  ~  f  -  (1  +  2E4>;  ‘  (8) 

in  which  case  there  exists  no  singularity  at  all  except  for  E  =  0, 
r  =  0.  Note  that  (8)  is  n^t  *  non-legitimate  solution.  A  special 
case  E  =  1  of  (8)  could  have  been  obtained  from  (2)  and  (3)  by  set¬ 
ting  v  +  X  =  0.  This  condition  indeed  yields  E  -  ±  1  hence  f  = 

±2 

(1  +  2$)  one  member  of  which  is  exactly  the  "standard"  line  element. 


The  other  member  is  just  as  un-objectionable  and  yields  the  three  ex¬ 
perimental  tests  correctly .  However,  it  has  no  singularity  and  no 
black  hoi el 

iv)  Consider  now  the  more  interesting  possibility ^ 

eV  =  f  =  1  -  2i)  +  2a<f>2  In  <}>  (9) 

where  a  is  an  arbitrary  number.  This  line  element  satisfies  the  New¬ 
tonian  limit  and  the  boundary  conditions  at  infinity  and,  by  the  same 
general  methods  of  calculation  mentioned  above,  yields  correct  values 
for  the  first  order  tests  of  red  shift  and  light  deflection.  For  the 
second  order  test,  namely,  the  perihelion  advance,  we  have 

=  (1  +  (10) 
h*  3 

which  is  a -dependent.  Does  this  mean  that  the  method  of  calculation 
is  not  as  general  as  it  was  supposed  to  be,  or  does  it  mean  that 
general  relativity  is  an  adjustable-parameter  theory  in  disguise? 
Furthermore ,  assuming  that  (9)  were  the  only  line  element  we  knew 
(other  being  assumed  not  yet  discovered) ,  how  does  one  calculate  a 
black  hole  behavior  free  from  ambiguities?  In  view  of  all  the  above 
does  the  usual  "null  surface"  characterization **  of  a  black  hole  have 
an  invariant  meaning? 


v)  Let  us  construct  the  Kretchmann  invariant  for  (4), 


I  =  ^Va\vaB  =  (12/m1*)  (1  -  f;G  (11) 

This  appears  to  assign  a  special  significance  to  the  "standard"  line 
element  as  I  =  12m2/r G  is  a  simple  power  of  r  only  in  that  case.  But 
if  this  is  of  any  invariant  significance  why  does  (11)  not  possess  a 
singularity  at  r  =  2m  as  an  invariant  characterization  of  a  black  liole? 
Note  also  that  it  is  not  possible  to  say  that  the  new  singularities  we 
presented  in  item  i)  are  merely  coordinate  singularities,  hence  they 
should  be  forgotten.  Because  in  view  of  all  that  is  shorn  we  can  equa¬ 
lly  well  argue  that  the  usual  singularity  r  =  2m  is  itself  a  coordinate 
singularity  hence  it  is  to  be  forgotten.  The  conclusion  seems  to  be 
that  despite  the  voluminous  literature  build  up  around  the  mysterious 
concept  of  a  black  hole,  Einstein' s  theory  does  not  really  predict  a 
unique  behavior  for  such  an  object. 

These  and  other  examples  that  one  can  easily  construct  [for  in- 

P  P 

stance  e°  =  f  =  1  -  2$  II  (1-B  ,<f>) ,  1  -  f  =  2$  II  (1  -  B  .<j>)  ,  B  .>  Oj 

j=l  3  j=l  3  3 

on  the  basis  of  the  general  solution  (4)  seem  to  lead  to  a  necessity 

of  examining  anew  the  status  of  the  problem  of  coordinates  in  general 
relativity .  5  One  could  say  that  due  to  the  coordinate  invariance  any  cal¬ 
culation  carried  out  in  the  standard  line  element  can  be  taken  over 
into  others  (by  transformation)  in  the  sense  of  redefining  the  inte¬ 
gration  variables  and  their  limits.  But  the  problem  is  that  we  do 
not  know  the  meaning  of  r  in  any  of  these  coordinates  (including  the 
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standard  coordinates)  and  we  do  not  know  the  integrations  limits. 

When  we  interpret  r  as  the  radius  in  the  standard  line  element  we 
assign  a  number  to  the  distance.  If  we  interpret  r  as  the  distance 
in  another  solution  we  should  assign  a  different  number  for  the  dis¬ 
tance  so  as  to  find  the  end  lesults  of  the  calculations  the  same. 

That  much  is  clear  and  trivial.  If  we  say  that  the  correct  assign¬ 
ment  of  a  number  is  possible  only  in  the  standard  solution  then 
either  we  are  attributing  an  ad_  hoc  privilege  to  the  standard  solu¬ 
tion  or  we  must  have  a  reason  for  the  correctness  of  this  assignment 
which  should  then  be  gcneralizable  to  other  solutions.  In  special 
relativity ,  from  which  the  general  theory  is  an  extension,  such  a 
problem  would  be  resolved  by  the  constancy  of  the  velocity  of  light, 

| dr |  =  c|dt|,  because  we  then  have  an  independent  way  of  measuring 
r  by  light  signals  in  "every"  Lorentz  line-element.  The  Lorentz 
line  element,  however,  is  not  generally  covariant  (for  example  we 
do  not  use  such  otherwise  covariant  substitutions  x'  =  x  +  Vt,  t’  =  t) 
but  covariant  only  under  such  substitutions  which  leave  the  velocity 
of  light  invariant.  It  seems  to  follow  that  a  similar  imposition  on 
(4)  must  be  made  to  identify  physically  the  meaning  of  r  by  an  inde¬ 
pendent  physical  means,  at  least  in  one  system  of  coordinates .  The 
general  theory  of  relativity  does  not  seem  to  satisfy  this  requirement 
although  everyone  believed  that  it  does.5 

The  conclusion  seems  to  be  that,  in  the  past,  special  non-covariant , 
choices  such  as  ]i  =  0,  V  -  \  -  0,  d^(V-g  g  )  =  0  tended  to  yield  at 


. . —  . .  . . . -i-r. 


least  a  plausible  way  of  imposing  the  extra  requirement  since  they 
eliminated  the  more  troublesome  choices  of  the  type  (9).  Similarly 

n 

the  usual  PPN  method'  which  relies  only  on  power  series  expandable 
metrics  imposes  an  extra  constraint  and  disregards  (9)  and  other  non¬ 
power  series  expandable  solutions.  Under  any  of  the  above  conditions 
I including  v  +  A  =  0  where  | £ |  =  1  or  ever  the  more  general  case 
| K |  <  2.5  x  205],  the  metric  (4)  yields  the  three  experimental  tests 
to  the  accuracy  of  the  present  observations  but  the  other  possible 
choices  are  still  too  large  in  number,  physically  non-equivalent  to 
each  other,  and,  do  not  answer  unambiguously  the  original  question 
we  posed:  Does  general  relativity  unambiguously  predict  the  existence 
of  black  holes? 
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APPENDIX  III 


PLANE  GRAVITATIONAL  WAVES 


a)  Plane  Waves  in  Einstein’s  Theory 

It  is  well  known  that  the  linearized  (first  order)  Einstein  theory 
has  plane  gravitational  waves  of  the  form1 

®00  =  ^  *  ^33  1 

gn  =  -1  +  4Q  g22  =  -1  -  4Q,  g1?  =  R  (2) 

where  Q  and  R  are  arbitrary  functions  of  t  -  z 

Q  -  Q(t-z)  ,  R  “  R(t  -  z)  (3) 

satisfying  D2Q  =  0,  □  2R  *  0.  It  must  be  stressed  at  once  that  for 

the  above  metric  C]2  =  D0^  reduces  to  □  2  ■  (82/8t2  -  82/3x2)  hence 

a 

Q  and  R  are  solutions  of  the  d'Alembert  equation  of  the  curved  space 
in  question.  Not?  further  that  forms  such  as  Q  =  t  z,  or  Q  = 

(t  -  z)a  are  not  acceptable  as  waves  since  they  do  not  satisfy  ap¬ 
propriate  wave  conditions.  The  solution  (2)  is,  as  is  well  known,  the 
basis  of  the  usual  calculations  of  the  energy  radiation2  and  the  ela¬ 
borate  experimental  investigations  to  detect  such  radiation.3 

In  the  course  of  our  investigations  it  became  clear,  however,  that 
the  solution  (2)  and,  together  with  it,  the  usual  belief  that  Einstein's 
theory  would  predict  gravitational  radiation  of  energy  gets  into  a 
peculiar  difficulty  in  second  order.  To  see  the  nature  of  the  diffi¬ 
culty  we  let 


-1  +  4XQ  +  X2U 


(4) 


g22  "  _1  +  4XQ  +  X2V 

gi2  -  4XR  +  X2W 


(5) 

(6) 


where  X  is  an  expansion  parameter  and  U,  V,  W  are  quantities  that 
are  of  second  order.  Substituting  into  the  Einsteinian  equations 


(7) 


one  finds  that  the  second  order  quantities  U  and  V  must  satisfy  the 
equation 


(U  +  V)"  +  32(QQ"  +  RR")  +  16(Q’2  +  R’2)  =  0  (8) 

where  prime  (')  implies  derivative  with  respect  to  t  or  z  or  £  =  t  -  z. 
(note  that  □  2f  =  0  is  not  equivalent  to  f"  =  0.)  Interestingly  W 
does  not  appear  in  the  equation  (8)  showing  that  the  non-linear  term 
in  g^  is  of  third  order  or  higher.  Being  a  differential  equation 
in  Q  and  R  equation  (8)  implies  that  U  +  V  is  expressible  as  a  func¬ 
tion  of  Q  and  R.  The  most  general  second  order  expression  is 

U  +  V  =  aQ2  +  0QR  +  yR2  (9) 


Substituting  into  (8)  one  sees  first  of  all  that  8=0.  Furthermore 

(crU6)QQ"  +  (cri-8)Q'2  +  (y+16)RR"  +(y+8)R'2  =  0  (10) 

Since  Q  and  R  are  linearly  independent  functions  it  is  sufficient 
to  analyze  the  case  of  Q  alone  (or  R  alone).  We  thus  set 

(a  +  16)QQ"  +  (a  +  8)Q'2  =  0  (11) 

A  trivial  solution  is  Q"  -  0,  Q  =  t  -  z  with  a  =  -8.  But  this  is 

clearly  unacceptable  because  Q  =  t  -  z  is  not  a  wave.  (If  Q  =  t  -  z 

were  acceptable  then  even  the  linearized  solution  (2)  is  not  unique, 
since  =  -1  +  aQ,  g22  =  -  1  +  bQ,  where  a  and  b  arbitrary,  would 
be  acceptable).  The  general  solution  of  (11;  is 

16+a 

Q  =  C  •  (t  -  z)24  +  2a  (12) 

This  Q  although  satisfies  the  d'Alembert  equation  □  2Q  =  0  is  never¬ 
theless  unacceptable  as  a  gravitational  wave  because  it  is  prescribed 
(Q  cannot  carry  information)  and  will  not  satisfy  wave  conditions  at 
infinity.  The  only  other  way  to  obtain  a  solution  is  to  set  in  (10) 
R  =  iQ  but  this  makes  R  and  Q  linearly  dependent  and  the  metr'f''  non- 

real.  The  conclusion  seems  to  be  that  in  Einstein's  theory  the 

linearized  solution  (2)  does  not  extend  to  second  order  consistently 
unless  Q  •=  R  =  0,  that  is  the  flat  metric. 
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b)  Plane  Waves  in  the  New  Theory 


In  the  new  theory  this  same  problem  goes  as  follows:  The  metric 


g  ■  n  e 
yv 


2($  -  2<j>) 


R  -Q 


where  <#>  =  =  0,n2$  *=  0.  Expansion  of  (13)  yields 


gu  «  -  cosh  4ft  +  ^  sinh  4ft 


g22  =  -  cosh  4ft  -  ^  sinh  4ft 


g12  “  o'  sinh 


800  *  "  g33 


where  ft  =  (R2  +  Q2)  .  This  metric  is  proven  to  be  an  exact  solution 

of  the  field  equations  of  the  new  theory.  The  wave  equation  is  indeed 
satisfied,  as  stated,  namely,  without  becoming  prescribed  and  the 
stress-energy  tensor  is  the  canonical  stress-energy  tensor  of  the 
Lagrangian 


(19) 


h  ■  -3X  3X*«  +  2  3**  V 


Since  $  ■  0  this  yields 


tV  =  -2  3  ♦“  3V$6  +  6V3X^6  3.$“ 
y  yT6  Ta  y  a  AS 


(20) 


For  R  “  0  the  solution  reduces  to  the  form 


-4Q  4Q 

*11  “  _e  »  &22  "  "e 


(21) 


whereas  for  Q  =  0  to  the  form 


gH  =  g22  =  "  cosh  4R>  g12  =  ±  sinh  4R 


(22) 


We  conclude  that  in  the  new  theory  the  linearized  solution  (2)  extends 

v 

consistently  to  second  order.  The  t^  then  provides  the  physical  basis 
of  energy  radiation  and  it  can  be  shown  that  the  amount  radiated  is5 


-  |£  -  ii-  D  TT 
3t  5C5  WV 


(23) 


where  is  the  quadrupole  moment  of  the  matter  distribution.  The 
implication  of  this  Appendix  seems  to  be  that  a  corresponding  calcula¬ 
tion  in  Einstein’s  theory  cannot  be  carried  out  consistently  because 
the  first  order  plane  wave  solutions  on  which  the  computation  depends 
do  not  carry  over  to  second  order.  The  second  order  consistency  is 
on  the  other  hand  absolutely  necessary  because  equations  themselves 
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v 

are  non-linear  and  the  stress-energy  tensor,  t^,  is  a  second  order 
quantity.  Consequently  if  experimentalists  discover  the  existence 
of  the  gravitational  waves  this  could  be  regarded  as  a  confirmation 
of  the  new  theory. 

Interestingly  Einstein  himself  seems  to  have  been  aware  of  the 
non  existence  of  gravitational  radiation  in  his  theory.  In  an  un¬ 
dated  letter  (1936-37)  to  Max  Bom  he  states:6  "Together  with  a 
young  collaborator,  I  arrived  at  the  interesting  result  that  the 
gravitational  waves  do  not  exist,  though  they  had  been  assumed  a 
certainty  to  the  first  approximation.  This  shows  that  the  nonlinear 
general  relativistic  field  equations  can  tell  us  more  or,  rather, 
limit  us  more  than  we  have  believed  up  to  now."  Unfortunately  later 
investigators  do  not  seem  to  have  paid  attention  to  Einstein's  con¬ 
clusion.  They  have  indeed  produced  exact  solutions  to  the  equations 
with  t  -  z  dependence  but  they  did  not  realize  that  these  solutions 
are  prescribed.  A  familiar  example  they  cite  is  of  the  form7 


4K  -  4Q 


4K  +  4Q 


4K 

where  K  is  a  second  order  quantity  (L  =  e  is  unity  in  first  order) . 
Expanding  as  in  (1)  (2)  and  (3)  one  finds  (K2  is  fourth  order  hence 
neglected) 


U  -  V  -  -8Q2  -  4K 


(25) 


Lotting  K  ■  eQ2  one  finds  a  =  -16  -  8e,  hence  the  problem  is  reduced 
to  the  previous  case  of  (12).  In  particular  if  e  *  -1,  a  *  -8,  hence 


Q  *  C  *(t  -  z).  This  case  is  not  only  prescribed  but  also  turns  out 
to  be  flat,  namely,  “  0. 


Another  example  cited  is  the  metric 


8yv  “  nyv  +  V  3vS  l(x2  "  y2>G(5)  +  H^>]  (26) 

where  G(£)  and  H(£)  are  arbitrary  functions  of  (  *  t  -  z.  The  problem 
with  this  solution  is  that  when  explicitly  written  out  it  yields 


g00  =  1  +  (x2  -  y2)  G  +  xy  H  =  1  +  A  (27) 

8q3  ■  -A  ,  g33  -  -1  +  A  (28) 

namely,  it  depends  only  on  a  single  function  A,  hence  gQ0,  gQ3,  g33 
are  not  linearly  independent.  To  represent  a  spin  2  graviton  propa¬ 
gating  with  the  velocity  of  light  one  however  needs  at  least  two 
linearly  independent  field  variables.  It  will  not  do  to  say  that  G 
and  H  are  linearly  independent  because  they  do  not  contribute  to  g^v 
in  a  linearly  independent  way.  At  the  present  time  we  are  not  aware 
of  the  existence  of  a  solution  of  Einstein’s  theory  which  can  be  con- 

t 

sidered  physically  as  a  gravitational  wave. 
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where  p2  =  A2  +  B2  +  C2. 
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(8) 


APPENDIX  V 


GENERAL  STATIC  CURVED  SOLUTIONS  OF  EINSTEIN'S  EQUATIONS  FOR 
SPHERICALLY  SYMMETRIC  METRICS* 


Richard  Pavelle 

Perception  Technology  Corporation,  Winchester,  Massachusetts  01890 


For  a  general,  static,  spherically  sym¬ 
metric  coordinate  system  we  find  the  general 
curved  solutions  of  the  Einstein  equations 
Gg  =  0  and  =  XS°.  trom  these  general  so¬ 
lutions  the  well  known  standard  solutions  are 
easily  generated. 


*This  research  is  supported  by  the  Advanced  Research  Projects  Agency 
of  the  Department  of  Defense  Under  Contract  No.  DAHC15  73  C  0369. 
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INTRODUCTION 

In  a  static  spherically  symmetric  coordinate  system  there  are 
a  number  of  well  known  particular  solutions  of  Einstein's  equations 
corresponding  to  specific  choices  of  the  energy -momentum  tensor. 
However,  the  problem  of  finding  the  general  solution  of  the  Einstein 
equations  corresponding  to  a  particular  matter  distribution  does 
not  seem  to  have  been  considered.  It  could  be  argued  that  in  view 
of  Birkhoff's  theorem1  this  problem  would  be  of  marginal  interest. 
This  is  not  the  case,  however.  On  the  one  hand  the  various  proofs 
of  the  Birkhoff  theorem  do  not  specify  a  method  for  constructing 
the  most  general  form  of  the  metric.  Secondly,  it  is  of  mathematical 
interest  to  find  general  solutions  for  any  system  of  non-linear  dif¬ 
ferential  equations  as  an  existence  proof  might  miss  some  parti¬ 
cular  solutions.  Finally,  although  the  general  solution  must  itself 
be  a  transform  of  some  particular  solution2,  there  is  no  a  priori 
guarantee  that  all  solutions,  for  a  given  matter  distribution  will, 
give  the  identical  predictions  of  physical  phenomena.  This  final 
point  is  discussed  elsewhere3  where  we  show  that  certain  transforms 
of  the  Schwarzchild  solution  do  not  seem  to  lead,  unambigously ,  to 
the  same  physical  predictions. 

Below,  we  solve  the  Einstein  equations  for  a  spherically  sym¬ 
metric  coordinate  system.  We  consider  both  the  vacuum  equations 
and  the  equations  with  the  cosmical  constant.  From  these  general 
solutions  certain  well  known  solutions  are  easily  generated. 

All  mathematical  statements  in  this  paper  have  been  verified 
and  intermediate  steps  simplified  using  the  symbolic  manipulation 
system,  MACSYMA,  developed  by  the  Math  Lab  group  at  Project  MAC, 
M.I.T. 


MaHUiMfc. 


wish  to  solve. 


We  approach  this  system  of  equations  by  noticing  that  G 


1 

1 


=  0 


may  be  solved  for  A  as 


A 


F' (DF'  +  2D"?) 
4DF 


(1.5) 


We  substitute  this  into  (1.3)  and  (1.4)  (set  to  zero)  to  find 


g3  =  D/F(2D/FF//  -  3D/(F/)2  -  2D/ /FF/)  _ 
3  (F/)2(DF/  +  2D'F)2 


(1.6) 


and 


g4  =  _  2D(2D/FF//  -  3D'(F')2  -  2D7 'FF')  _ 

4  (F ' ) (DF '  +  2D'F)2 

It  is  clear  that  if  D  4  0  and  F  ^  0  the  most  general  way  we  may 
satisfy  (1.6)  and  (1.7)  simultaneously  is  by  requiring 


2D,FF/ '  -  3D ' (F ' ) 2  -  2D/,FF/  =  0 


(1.8) 


This  differential  equation  is  easily  integrated  to  yield 


D  =  K 


1 


where  Kj  and  K2  are  constants. 


(1.9) 


. 


In  (1.6)  and  (1.7)  we  note  that  F'  ^  0  and  DF'  +  2D'F  ^  0.  These 


"restrictions"  on  the  functions  are,  from  (1.5),  nothing  more  than 
the  requirement  A  ^  0  which  we  stipulated  above. 


We  now  substitute  (1.9)  into  (1.5)  to  obtain 


K^F')2 


4F*  (K.J  -  K2F'1/2) 


(1.10) 


From  (1.9)  and  (1.10)  in  addition  to  a  trivial  coordinate  transforma¬ 
tion  we  then  have5 

THEOREM:  The  metric  tensor  defined  by  the  line  element 


ds2  =  - 


(FT 


4F(1  -  KF_1/2) 


dr2  -  F(d02  +  sin20d<|>2) 


+  (1  -  KF_1/2)dt2  ,  (1.11) 


satisfies  the  Einstein  vacuum  equations  identically  where  _F  is^  an 
arbitrary  function  of  r_  of_  class  C2. 

We  may  transform  (1.11)  into  a  more  intuitive  form  by  letting 
F  -*  B2r2  where  B  =  B(r)  is  arbitrary.  With  this  new  definition  it 
is  easily  seen  that 

2 

ds2  =  -  ■  dr2  -  B2(r2d02  +  r2sin20d$2)+(l  -  |-)dt2  . 

Br 


(1.12) 


The  non-isotropic  Schwarzchild  solution6  may  now  be  obtained  by  set¬ 
ting  B  =  1  to  find 


F 


ds2  - - dr2  ‘  r2(d02  +  sin26d$2)  +  (1  -  -)dt2  .  (1.13) 

1  -  -  r 


From  (1.11)  we  may  also  obtain  Fock's  solution7.  By  choosing 
(r  +  K/2)2  we  find 


ds' 


♦! 


K 

r  "  2 


dr2  - 


(r  +  |)2(d02 


+  sin20d<}>2)  + 


K 

r  -  2 


dt2.  (1.14) 


From  (1.11)  or  (1.12)  we  may  find  the  isotropic  spherically 
symmetric  line  element  by  setting 


(F')2 


4F(1  -  KF_1/2)  r2 


(1.15) 


This  differential  equation  may  easily  be  solved  for  F  and  we  find 

(1.16) 


+f/r2 


4K2 


where  Kj  is  an  integration  constant.  Thus,  the  isotropic  line  element 


takes  the  form 


(  i-JL)2 

da2  »  -  (1  +  7-  )  (dr2  +  r2d02  +  r2sin20d<|>2)  + -  dt2  (1.17) 


In  accordance  with  the  well  known  result.8 


2.  EQUATIONS  WITH  THE  COSMICAL  TERM 

For  the  metric  (1.1)  it  follows  from  (1.2),  (1.3)  and  (1.4) 
that  the  component  form  of  the  Einstein  equation 


G 


i 

j 


(2.1) 


where  X  is  the  "cosnical  constant"  is 


G1  =  D(4AF  -  (F7)2)  -  2FF7P7 
1  4AF2D 


(2.2) 


G2  =  g3  =  D2tA(2FF/'~(F/)2)~A/FF/]+D^A(D/FF,+2p/,F2)~A/D,F21~AF2(D/)2 
2  3  4A2F2D2 


X  (2-3) 


,4  4A2F2  +  2A7FF7  +  A((F7)2  -  4FF77) 

a - -  _ 

4A2F2 


(2.4) 


As  we  did  for  the  more  simple  case  we  solve  (2.2)  for  A  to  find 


F7  .  DF7  +  2D7F  - 
4FD  {  FX  -  1  } 


(2.5) 


ho  ho 


In  both  cases  we  then 


We  substitute  this  into  (2.3)  and  (2.',.  . 
subtract  (the  right  hand  side)  from  the  left  to  obtain  differen¬ 
tial  equations  which  turn  out  to  be  similar  to  each  other.  This 
procedure  gives 


and 


It  is  clear  that  to  solve  (2.6)  and  (2.7)  we  need  to  find  the 
relation  between  F  and  D  which  causes  the  expression  in  curly  brac¬ 
kets  to  vanish.  For  A  =  0  we  know  that  the  relation  we  seek  must 
reduce  to 

D  =  Kj  -  K2F'1/2  (2.8) 

hence  we  now  seek  a  function  N(F)  such  that 

D  =  Xj  -  K2F'1/2  +  AN(F)  (2.9) 


satisfies  f2.6)  and  (2.7). 
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Substitution  of  (2.9)  into  (2.6)  results,  after  some  laborious 
calculation,  is  the  following  differential  equations  for  N: 


N  -  N'F  *  0 


3N'  -  Kj  -  0 


(2.10) 


These  are  easily  integrated  to  find 


N 


K1F 


(2.11) 


Thus,  (2.9)  becomes 


D  =  Kx(l  -  f-)  -  K2F"1/2 


(2.12) 


With  this  expression  for  b  we  find 


A  = 


K^F')2 


4(Kj  -  K2F_1/2-  y  KjF)F 


(2.13) 


From  (2.11)  and  (2.14)  and  a  trivial  coordinate  transforma¬ 
tion  to  absorb  an  integration  constant  we  have5 

THEOREM:  The  metric  Refined  by 


ds2 


(F')‘ 


—\/2  \ 
4*F(1-KF  '  -  ±  F) 


dr2  -  F(d62  +  sin20d$2)+(l  - 


KF"1/2-  -|F)dt2 


(2.14) 
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where  C  is  an  integration  constant.  This  integral  does  not  appear 
to  be  expressible  in  closed  form.  However,  one  can  obtain  expres¬ 
sion  for  particular  cases  of  interest.  By  setting  K  -  0  in  (2.18) 
we  may  perform  the  integration  and  solve  for  F  to  find 


12  CR2 
(XR2  +  C)2 


(2.19) 


'Then 


ds2  - -2--  (dr2  +  r2d02  +  r2sin20d<f>2)  +  (1 - — -)dt2  (2.20) 

(XR2  +  C)Z  (XR2  +  C)2 

Here  we  have  the  isotropic  static  form  of  the  solution  of  (2.1).  Other 
solutions  of  possible  cosmological  interest  could  be  found  by 
performing  a  series  expansion  on  the  left  hand  side  of  (2.18). 


We  have  not  considered  such  cases  at  the  present  time. 
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APPENDIX  VI 


SOLUTIONS  WITH  PLANE  SYMMETRY  IN  EINSTEIN'S  THEORY 


Consider  the  Cartesian-like  coordinate  system  defined  by 


ds2  =  -  F(dx2  +  dy2)  -  Hdz2  +  K  dt2  (1.1) 


where  F,  H  and  K  are  function  of  z.  We  wish  to  find  the  general 
relations  between  these  functionals  so  that  (1.1)  satisfies  Einstein's 
Vacuum  Equation: 

For  (1.1)  we  find  the  components  of  the  Einstein  tensor  (set  to 
zero)  to  be 


G1  =  G2  =  -  {K(H(2F2K"  +  FF ' K '  ]  -  F2H'K')  -  F2H(K')2 

+  K2(H(2FF"  -  (F')2)  -  FF'H')}  /  4F2H2K2  =  0  (1.2) 


G3  =  _  F'  (2FK'  +  jgll  =  o 

3  4F2HK 


'1.3) 


G4  =  -  H(F[2F2F"+F(F’)2]  +  F2(2FF"-(F')2)  -  F2(F')2)-2r3F'H' 
4  ’  "  4F4H2 


(1.4) 


We  may  satisfy  (1.3)  by  requiring 


2FK'  +  KF'  =  0 


(1.5) 


This  differential  equation  may  be  integrated  and  we  find 


F  -  Cj/K2  (1.6) 

where  Cj  is  a  constant.  This  relation  implies  that  (1.2)  and 
(1.4)  become 

G1  =  G2  =  H(2KK"  -  5(K')2)  -  KH'K'  _  Q  (1  ?) 

1  2  ’  4H2K2 


H(2KK' '  -  5(K')2)  -  KH'K' 


(1.8) 


Hence,  by  solving  (1.7)  or  'a. 8)  we  will  have  a  general  sjiution 
of  the  Einstein  Vacuum  Equation.  These  are  integrablc  and  we  find 
that 

(K')2 

It  -  C.  -  (1.9) 

2  K5 

satisfies  (1.8)  and  (1.8).  Thus,  the  metric  tensor  components 
defined  by 
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On  the  otner  hand,  we  may  satisfy  (1.3)  by  choosing  F'  *  0. 
With  this  choice  it  is  found  that 


K"  +  -K  H/K/  -  H(K'): 


(1.11) 


The  solution  of  (1.11)  is 


(1.12) 


CjOC')2 


(1.13) 


We  thus  find 


ds2  =  -  (dx2  +  dy2)  - 


(K')2 


dt2  +  K  dt2 


(1.14) 


also  satisfies  G7  =  0. 

J 

However,  it  is  not  difficult  to  prove  that  for  the  metric  (1.14), 
all  components  of  the  Riemann  Christoffel  Tensor  vanish  identically: 
Hence,  the  metric  (1.34)  is  flat.  It  can  also  be  shown  that  (1.14)  is 
the  only  flat  metric  belonging  to  the  general  class  (1.1). 


From  (1.10)  we  may  seek  the  general  isotropic  line  element  by 


setting 


cl  c2(k')2 


(1.15) 


Then, 


K'  -  C  K3/2 


or 

K  »  (Cx  +  D)2 


Hence 


ds2  *  -  (Cx  +  D)  \dx2  +  dy2  +  dz2  ) 
+  (Cx  +  D)2  dt2 

is  the  isotropic  line  element. 


(1.16) 


(1.17) 


(1.18) 


. . . 
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FIELD  THEORY  OF  GRAVITATION * 


HUseyin  Yilmaz 

Perception  Technology  Corporation ,  Winchester ,  Massachusetts  01890 


SUMMARY 


A  field  theory  of  gravitation  where  the  stress- 
energy  tensor  of  gravity  contributes  to  geometric 
curvatures  is  described.  The  field  is  locally 
Lorentz  covariant,  namely,  in  an  infinitesimally 
small  laboratory  (not  necessarily  freely  falling) 
the  equations  are  formally  the  same  as  they  would 
be  in  special  relativity.  The  theory  reproduces 
all  known  gravitational  effects  correctly  and  has 
experimental  consequences  differing  from  those  of 
the  Einsteinian  theory. 


1.  PRELIMINARY  CONSIDERATIONS 


We  propose  to  construct  a  field  theory  of  gravitation  where 
special  relativity  is  satisfied  as  a  correspondence  limit.  To  make 
the  essence  of  this  correspondence  explicit  we  introduce  two  postu¬ 
lates:  a)  Matter  tensor  au^uv  generates  a  gravitational  field 
which  is  a  local  second  order  tensor,  b)  Local  velocity  of  the  gravi¬ 
tational  field  is  equal  to  c.  We  express  the  content  of  these  two 
statements  as 

uyuV  (1.1) 


9  V 

=  47 ra 
y 


3 

v  y 


0 


(1.2) 


It  is  evident  that  the  two  postulates  are  generalizations  of  their 

special  relativistic  counterparts  in  the  context  of  gravity.  The  first 

postulate  extends  one  of  the  most  important  consequences  of  special  re- 

v 

lativitj,  namely ,  the  matter  density,  a,  becomes  a  tensor  a u  u  .  The 

V  V  o  ^ 

second  postulate  requires  3^$^  =  0  otherwise  a  A3  3a<^  could  be  added 
to  the  left  hand  side  of  (1.1)  which  alters  the  vacuum  velocity  of 
propagation  for  some  components.  Locality  further  implies  that  the 
absolute  values  of  the  fields  are  unobservable  and  that  only  the  dif¬ 


ferences  $  (x)  -*■  (x)  -  $(x' )  enter  into  the  observable  phenomena.  Note 

that  (1.2)  may  further  be  supplemented  by  3  3^$^  =  0  because  when  one 
contracts  (1.1)  the  additive  terms  must  still  be  absent  from  the  propa¬ 


gation  equation  of  the  trace,  $  =  The  gravitational  field  so  con- 

^  v 

ceived  must,  of  course,  have  a  stress-energy  tensor  t^  which  must  add 
to  the  matter  tensor  as 


v  v  v 

T  =  o  u  u  +  t  /4i r  (1.3) 

y  y  y 


and  the  equations  of  motion  of  test  particles  must  be  obtained  as 


d2  x 

o - —  =  3  (t>/4v)  (1.4) 

ds2  V  y 


both  of  which  are  prerequisites  of  having  a  local  field  theory  in  cor¬ 
respondence  with  special  relativity. 


As  can  be  seen  these  equations  would  be  contradictory  within  the 


flat  space-time  of  special  relativity  because  the  usual  conservation 
v 

law  3^2^  =  0  would,  on  account  of  (1.2),  lead  to  separate  conservation 
of  matter  and  field  counterparts.  This  would  mean  that  there  can  be 


no  interaction  between  matter  and  field,  hence  (1.4)  would  be  meaning¬ 


less.  This  difficulty  arises  because  in  flat  space  3^  commutes  with 
□2.  Clearly  in  adopting  the  above  equations  we  have  implicitly  assumed 


a  more  general  space-time  where  we  require 


avn2  i<n2av  n.5) 

The  latter  is  realizable  in  a  curved  geometry. 


2.  DEVELOPMENT  OF  THE  THEORY 

In  the  curved  space-time  the  total  stress-energy  tensor  will  be 
identified  as  the  divergence- free  tensor 1 

V  1  .V  V  V 

R  -  =;  5  R  -  8u(du  u  +  t  /4t\)  (2.1) 

U  2  u  U  V 

Upon  contraction  we  shall  write  R  as  a  scalar  relation 

R  -  -8m  +  4L.  +  43°J  (2.2) 

t  a 

where  =  -t/2  is  the  field  Lagrangian  and  9aJa  a  divergence.  Our 
aim  is  to  determine  and  9°^  by  comparing  (2.2)  with  the  geometric 
expression  of  R,  namely, 

R  =  -  c9ln  Agr  -  |  }  9agryv  -  (J-g)~l  gvv  )  (2.3) 

Contracting  (1.1)  we  first  obtain  the  equation  of  $ 

□  2<fi  =  4m  (2.4) 

which  implies  that  in  (2.3)  we  must  have 

J-g  =  e ^  (2.5) 

Then  the  last  term  of  (2.3)  (which  is  the  only  one  that  can  contain  the 
Lorentz  condition)  and  (1.2)  imply  that  /-g r  gr^V  is  independent  of  <j>. 


In  view  of  (2.5)  this  requirement  implies,  in  tarn,  that  gr 
the  form 


is  of 


9 


n 


(2.6) 


where  $  is  the  symmetric  field  tensor  By  expanding  this  metric 

in  first  and  second  order  we  have 

g  =  n  (l  +  2$  +  242)  -4 n  +  an  a<J>a<f>6  (2.7) 

yv  yv  T  r  yXYv  Y  yX  v  ct0  y  v 

hence,  to  the  required  order,  the  last  two  terms  of  (2.3)  are 

- 1  <;vv  -  *  sv«  <2-s> 

-  (S-9)~h*v(S-g  gyV;  =  -43y3v^  -  a3g*V<^  (2.8') 

We  therefore  find  L.  and  3aj  as 

f  a 

Lf  ■  3x*a  +  7  3**V  l2-9> 

33Ja  =  -  3y3v  =  0  (2.10) 

Concerning  the  equations  of  motion  (1.4)  it  is  not  difficult  to 
see,  from  the  form  of  (2.9),  that  the  canonical  stress-energy  of 
will  yield  the  local  force ^ 


a 


d2x 

_ y 

ds 2 


8  oi 

-2au  u  3  +  o3  A 

ct  yY0  y 


(2.11) 


Upon  dropping  a  these  equations  are  found  to  be  the  geodesic  equations 
of  the  space  given  by  (2.6),  namely. 


d2x 

_ y 

ds 2 


1  ,  a  0  ,  a-,  £ 

=  —  3  gr  .u  u  ={.}uu 

2  y’o0  y0  a 


(2.12) 


Thus,  under  the  condition  that  is  an  inertial  space  the  presence 
of  gravity  leads  to  a  curving  of  that  space  as  in  (2.6).  Moreover, 
to  the  approximation  here  described,  the  gravitational  equations  of 
motion  are  the  same  as  the  geodesic  equations  of  the  curved  spare  so 
generated. 


3.  APPLICATIONS 


In  the  usual  solar  system  applications  the  source  is  static  so 

V  0  0  0 

that  -*•  Since  in  this  case  <j>g  =  4>,  3  <f >  =  0  the  problem  is  very 

simple.  From  (2.6)  we  have 


9 


00 


(3.1) 


where  <f>  =  M/r.  The  form  of  gQ Q  and  g^  show2 that  the  red  shift  and 
light  bending  are  correctly  predicted.  As  for  the  perihelion  advance 
we  transform  (3.1)  into  polar  coordinates.  Since  in  this  case  <f>  is  a 
function  of  r  only,  the  qeodesic  equations  of  motion  (2.11)  yield  the 
constants  of  the  motion,  un  -  k,  u  =  q,  u  =  h.  Choosing  0  =  ir/2, 

U  O 

uQ  =  0  one  has 


u 


0 


(3.2) 


where  k  -  1.  The  line  element  now  is 

ds2  =  e  2^dt2  -  e2^ (dr2  +  r2dy2)  (3.3) 


Upon  introducing  (3.4) 


into  (3.5)  and  eliminating  dt  and  ds  one  has 
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y 

\ 


Setting  u  =  1/r,  and  differentiating  one  finds 

U  +  u  »  f(u)  =  —  e2* (2e2*  -  1)  (3.5) 

h2 

which,  as  is  well-known,2  yields  the  observed  perihelion  advance, 

A  ft  -  vdf/du  ]  =  6tiM2/h2,  a  =  f(0)  =  M/h2 .  Likewise  a  calculation  of 

u-a 

the  radar-echo  delay  with  (3.1)  and  (3.2)  gives  the  observed  value 
(l  and  l'  are  the  orbital  radius  of  the  planet  and  the  earth,  re¬ 
spectively)  . 


At  =  4M 


f 

tin  4 


(3.6) 


to  within  the  experimental  accuracy. 

In  the  weak  quasistatic  limit  where  -  <j)  is  the  largest  compo¬ 
nents  (2.6)  yields 


g  -  n 
mV  yv 


(l  +  2$  +  2 *2;  - 


(3.7) 


A  comparison  with  the  usual  PPN  formalism  then  shows  that  all  the  pre¬ 
ferred-frame  parameters  are  zero  hence  no  preferred  frame  effect  is 
predicted.  'J'hs  theory  therefore  reproduces  correctly  all  the  experi¬ 
mental  results  concerning  earth  tides  and  to  the  absence  of  certain 
effects  which  would  otherwise  exist  in  a  preferred  frame  theory.  Thus 
at  first  sight  the  predictionsof  the  new  theory  seem  to  be  the  same  as 
those  of  Einstein' s  theory .  ^ 

There  are,  however,  situations  in  which  the  present  theory  differs 
drastically  from  Einstein's  theory.  For  example  the  vanishing  covariant 
divergence  of  (2.1)  yields 


(i/-g)  *3  (J-g  au  uV )  -  a 
v  y 


{  u  +  D  (t^ / 4^)  =  0 
y6  a  v  y 


(3.8) 


Since,  in  first  and  second  order,  the  last  term  is  the  same  as  in 
(1.4)  one  finds 


\(S-9  auyu^  =  0  (3.9; 

which  expresses  the  conservation  of  total  energy-momentum  of  matter 
in  the  presence  of  gravity.  It  is  well-known  that  Einstein's  theory 
yields  d^(S-g  Qu' )  =  0  which  is  the  conservation  of  total  rest-mass . 
Since  in  the  limit  of  special  relativity  (3.9)  reduces  correctly  to 
the  conservation  of  energy-momentum  of  matter,  we  are  able  to  accom¬ 
modate  the  increase  of  rest-m^cs  when  two  objects  coalesce  in  a  ra¬ 
diationless  collision  (the  so-called  mass-defect)  and  thereby  preserve 
m  =  E/c2.  The  argument  is  in  fact  valid  both  in  the  strong  and  in 
the  weak  field  limits  and  allows  e+  +  e~  -►  2y.4 

In  the  extremely  strong  field  limit,  near  Schwarzchild  radius, 

the  behavior  of  (3.1)  is  seen  to  be  drastically  different.  In  fact 

there  is  no  singularity  except  at  r  -*•  0  hence  no  black  hole  behavior 

4,5 

in  the  usual  sense  is  predicted.  ‘.’he  collapsed  object  will  here 
remain  as  naked  except  for  r  =  0.  In  the  case  of  neutron  stars  with 
radii  considerably  larger  than  a  =  GN/c2  the  static  behavior  is  not 
very  different  from  the  predictions  of  the  conventional  theory.  As 
for  the  highly  non-static  cases  the  theory  predicts  gravitational 
radiation  from  exploding  galaxies,  which  might  be  detectable . 6  In 
contrast  the  Einsteinian  theory  seems  to  be  unable  to  accommodate 
gravitational  radiation.  In  a  recent  calculation  with  the  help  of 
powerful  symbolic  manipulation  techniques  it  is  shown  that  the  con¬ 
ventional  linearized  plane  wave 

gu  -  -1  +  4Q  ,  g22  =  -1  -  4Q  ,  gl2  =  4R  (3.10) 

where  R  and  Q  are  functions  of  t  -  z  would  lead  to  contradictions  in 
second  order.  More  specifically  it  is  found  that  only  three  special 
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cases  admit  a  solution  a)  Q  =  R  =  constant  (flat  space) ,  b)  R  =  iQ 

(metric  is  not  real),  c)  R  =  0  (no  polarization)  g..  =  , 

4  (K-0) 

=  -e  where  K  is  second  order.  Substituting  into  R  =  G 

22  X  yv 

one  finds  Q  =  C*  (t-z)  .  ( Unacceptable  at  z  ■*  <*>;  Q  prescribed) . 


2  1  2 

.Tn  the  new  theory  it  is  seen  that  $„  =  -$,=  Oft-zJ,  (J>,  =  R(t-z) 

Z  X  * 

yields  the  metric  g^  -  -  cosh  4(1  +  (0/d)  sinh  4d,  g22  =  -  cosh  4(1  - 
(0/d)  sinh  4d,  g12  =■  (R/d)  sinh  4d,  gQ0  =  “  #33  =  y ields  a  solution 


to  the  field  equations  hence 


*  = 


/• 

! 

!  . 


1  • 

1 

V 


\ 

*  *\ 
•  • 

0  R 
R  -Q 


(3.11) 


is  a  plane  gravitational  wave.  Interestingly  this  leads  to  (3.10)  as 
its  first  order  approximation.  Likewise  it  can  be  shown  that  in  Ein¬ 
stein's  theory  a  first  order  solution,  which  is  necessary  in  view  of 
a  200  km/sec  motion  of  the  solar  system  relative  to  the  center  of  ga¬ 
laxy,  becomes  unsatisfactory  in  second  order.  Again  the  new  theory  is 
free  from  the  difficulty. 


4.  POSSIBLE  EXTENSIONS 


The  present  theory  can  be  completed  with  regard  to  mechanical  forces 

such  as  pressure,  and  forces  of  electromagnetic  and  nuclear  nature,  by 

following  a  well-known  procedure.  We  shall  not  repeat  this  process 

which  is  more  or  less  straightforward.  Instead,  we  shall  present  a 

more  speculative  possibility  which  resembles  the  Einsteinian  idea  of 

a  unified  field  theory  whereby  various  fields  are  introduced  into  the 

v 

framework  in  analogy  to  $  .  Let  (2.6)  be  written  as 


g  =  n 


2($  -  2$)  +  2z 


(4.1) 


where  e  is  an  antisymmetric  tensor.  We  define  a  tensor  f  such 

yv 

that  e  is  the  world-line  integral  of  f 

yv  3  yv 


yv 


iff  ds 
2  yv 


(4.2) 


Computing  R  it  can  be  shown  that  the  Lorentz  conditions  are  automa¬ 
tically  satisfied 


3P3  fV  =  0 
v  y 


If  we  define  aconserved  quantity,  j 


3  f  =  -  4tt7 
v  y  “'y 


3Pj  =  0 

y 


(4.3) 


(4.4) 

(4.5) 


we  have  the  first  set  of  Maxwell's  equations  plus  the  conservation  of 
charge.  Then  the  antisymmetry  of  £  yields  the  second  set 


3  f  +  3  £  +3  £ 

v  ya  a  vy  y  ov 


(4.6) 


It  can  further  be  shown  that,  by  a  suitable  definition  of  {  “},  the  sign 

y3 

of  the  electric  charge,  the  acceleration  caused  by  t,  becomes 


dzx 


us 


—  =  3  (tV/4^ i)  =  £  jV 

„2  v  y  yv^ 


(4.7) 


This  leads  to  a  Lorentz  force 


d2x 


m 


ds 


u  ,  v 
—  =  ef  u 
2  yv 


(4.8) 


where,  by  integration,  o  and  p  are  replaced  by  the  rest-mass  and  the 
charge.  We  can  see  that,  in  at  least  a  formal  sense,  equations  si¬ 
milar  to  those  of  the  electrodynamics  are  obtainable  from  such  a 
process  as  a  local  field  theory.  However,  the  t  of  the  non-symmetric 
part  does  not  reproduce  the  Maxwell  stress-tensor  exactly.  We  shall 
point  out  only  that,  by  a  suitable  introduction  of  magnetic  poles. 


hopes  are  recently  entertained  that  an  SU(3)  symmetry  for  heavy  par¬ 
ticles  might  be  obtainable.7  This  direction  may  conceivably  be  a 
remedy  to  the  above  di-ficulty  with  the  Maxwell  tensor.  The  problem 
is  under  further  considerations. 

Less  speculative  is,  perhaps,  the  possibility  of  a  pure  cosmical 
scalar,  of  Dicke-Brans  type.  Again  writing 


9 


yv 


e2(b  -  2$;  +  X* 

yv 


(4.9) 


where  X  is  a  small  admixture  parameter,  we  can  see  that  such  a  field 
can  only  yield  a  conformal  extension  of  the  present  theory.  Further¬ 
more  if  X4  >  0.02<J>  one  has  to  renormalize  sr0g  in  order  not  to  conflict 
with  the  red-shift  experiment.  A  pure  scalar  of  this  type  to  repre¬ 
sent  the  whole  gravitation  would  of  course  be  out  of  the  question  be¬ 
cause  it  would  not  produce  any  bending  for  light.  For  a  small  admix¬ 
ture  X  =  0.05  'V'  0.08  the  bending  of  light  and  radar  echo  delay  are 
correspondingly  reduced  and  the  theory  can  be  tested  by  experiment. 

The  essential  point  is  that  such  a  field  is  possible  as  a  conformal 
extension  of  the  present  theory.  An  extension  of  a  different  kind, 
namely,  the  possibility  of  a  rest-mass ,  y,  for  gravity  seems  to  be 
unrealistic  on  account  of  the  very  long  range  of  the  gravitational  in¬ 
teraction. 


5.  COVARIANCE 

v 

In  Section  2  we  have  seen  that  the  Lorentz  condition  3  <t  -  0 

v  y 

leads,  in  first  order,  to  the  expression,  gVV )  =  0,  which  is 

reminiscent  of  the  so-called  harmonicity  conditions.  Since  the  usual 
harmonicity  conditions  are  non-covariant  a  question  arises  as  to  the 
nature  of  the  coordinate  invariance  satisfied  by  the  new  theory.  We 
shall  now  show  that  although  they  reduce  to  each  other  in  the  trivial 


teiiffuitMma.-wt-mi.-Mri  ■■■«■<  nmiurr ii'rtriiiifi 
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case  of  cartesian  coordinates  our  conditions  and  the  usual  harmonicity 
concept  have,  in  general,  nothing  to  do  with  each  other. 

Let  a  curved  space  h  (x)  be  defined  over  all  curvilinear  coordi- 

yv 

nates  (x)  of  special  relativity  such  that 


h  (x)=t\  (x)g.  (x) 

yv  yX  3Xv 


(5.1) 


where  n  =  (1,1, 1,1) ,  n  =  (1,  r2,  r2sin2d,  1)  etc.  The  g^  (x)  are 
functions  of  space-time  variables  and  they  act  as  congruent  coeffi¬ 
cients  to  Upon  constructing  the  d'Alembertian  of  h^(x)  one 

finds 


n 2  =  gryX(V nJ  nXv3y;  +  (J-g)  ld^(S-g  gryv;3t 


(5.2) 


Since  the  first  term  of  this  expression  is  a  linear  combination  of 
special  relativistic  d ' Alembertians  the  local  propagation  velocity 
will  be  c  if  the  second  term  vanishes.  We  have8 


3 v^-sr  =  0 


(5.3) 


These  relations  are  not  coordinate  conditions .  On  the  contrary 
they  are  the  requirements  that  the  local  velocity  of  propagation  be 
c  as  in  special  relativity .  Subject  to  this  statement,  which  is  one 
of  our  principles ,  the  theory  is  covariant,  namely,  it  allows  all 
curvilinear  substitutions  permissible  by  the  local  universality  of  c. 
Note  further  that  the  conditions  (5.3)  do  not,  in  general,  imply  the 
harmonicity  of  the  coordinates  (for  example  in  spherical  coordinates 
corresponding  to  (3.1)  G2r  =  -  2e  v /r  ?  0) .  It  would,  therefore, 
be  quite  inappropriate  to  call  (5.3)  harmonicity  conditions.  We  sug¬ 
gest  that  they  be  called  the  "propagation  conditions"  since  they  have 
to  do  with  the  local  propagation  velocity ,  c.  It  is  clear  that  the 
theory  here  proposed  does  not  violate  the  principle  of  covariance . 
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The  coordinates  (2.6)  are  not  preferred  coordinates.  They  are  special 
only  in  the  sense  that  they  allow  a  local  special  relativistic  field 
theory  interpretation  of  gravity  and  a  local  special  relativistic 
philosophy  of  measurement. 

6.  EQUIVALENCE 

It  must  not  be  supposed  that  (2.6)  and  (2.12)  can  only  accommodate 
curved  spaces  and  their  fields  of  acceleration.  We  can  show,  with  the 
help  of  (5.3),  that  the  theory  accommodates  accelerated  frames  of  re¬ 
ference  in  an  otherwise  flat  space.  Let  n  =6  and  let  h  ( x )  be 

r  liv  yv  yv 

still  a  flat  space.  Since  R  ^  =  0,  everywhere ,  we  also  have  R  =  0, 
everywhere.  Thus  in  (2.2)  a  =  0,  Lf  =  0,  =  0.  Letting  the  acce¬ 

leration  to  be  uniform  in  the  z-direction  (2.12)  yields 

<j>°  =  -  yz  +  C  (6.1) 

3  o 

Then  the  propagation  conditions  (5.3)  are  satisfied  if  <t>3  =  Q0  =  4i. 

This  field  indeed  also  satisfies  the  other  conditions  such  as  0  =  0, 

R  =  0  etc.,  above,  hence 

-2  ill  -2ip 

?00  =  6  '  ”'32  =  e  '  =  ?22  =  "  1  {6'2) 

is  a  flat-space  solution  consistent  with  (2.6).  It  may  be  conjectured 
that  this  solution  represents  a  kinematically  accelerated  frame  of  re¬ 
ference  in  an  otherwise  flat  space  of  special  relativity . 

Note  that  if  the  conjecture  is  valid,  namely,  if  the  kinematical 
acceleration  does  not  curve  the  space,  as  gravity  would,  then  the  equi¬ 
valence  of  gravity  and  inertia  will  not  be  exact  but  restricted  to  the 
Newtonian  limit  v  0  alone.  In  fact  the  two  solutions  (3.1)  and  (6.2) 
yield,  respectively ,  the  geodesic  accelerations 


\ 


(6.3) 


2  v 

-it  =  3  $[i  +  —  ],  -it  *  3  ip[2  -  —  ] 

V  V  2  V  v  ' 

Cz  C“ 

For  v  0 ,  these  accelerations  are  equal  but  for  v  c  they  are  dif¬ 

ferent.  The  conjecture  seen  '  to  be  testable  by  an  orbiting  gyroscope 
placed  in  a  spherical  shell  Assuming  <J>  =  $  the  gyroscope  will  have, 
relative  to  the  shell,  a  differential  acceleration 

v2  +  v2  +  2V2 

-  De  =  3  *  — - j* - -  (6.4) 

v  v  6 

c 

The  gyroscope  could  conceivably  be  replaced  by  some  atoms  with  high 
rotational  energy  or  with  light  waves.  Since  for  a  good  gyroscope 
the  factor  is  of  the  order  of  !0-12  it  is  also  conceivable  that  a 
Dicke-Braginski  type  experiment  might  turn  out  to  be  feasible. 

The  outcome  of  such  an  experiment  would  be  of  great  importance 
to  see  if  the  theory  of  gravitation  has  to  incorporate  inertia  only 
in  the  limit  of  v  0  or  more  generally  for  all  velocities.  In  the 
latter  case  the  flat  space  solutions  of  the  kind  (6.2)  would  be  dis¬ 
carded  from  the  theory  by  saying  that  under  kinematic  accelerations 
space- time  becomes  curved.  On  the  other  hand  if  the  force  given  by 
(6.4)  exists  the  kinematic  accelerations  are  representable  by  coor¬ 
dinate  transformations  ■-  0)  whereas  gravitational  accelerations 

are  not.  In  fact  if  (6.4)  exists  the  equations  of  motion  (1.4)  will 

have  to  be  restated.  Thus  if  n  (x)  is  flat  but  nevertheless  non- 

y  v 

inertial  the  geodesic  equations  of  motion  (2.12)  will  contain  inertial 
forces  coming  from  as  well  as  the  gravitational  forces  arising 

from  g  ($) .  Thus  we  arrive  at  a  more  general  conception  of  space 

V 

V  WnM  9(h) vv  -  9vv(x.h  (6.5) 

where  <f>  =  0,  h ^  =  g  (x)  contains,  in  general,  arbitrarily  accelerated 
frames  such  as  rotation  and  accommodates  even  intrinsic! y  curved  spaces 


such  as  the  presence  of  a  cosmological  constant,  A.  The  presence  of 

gravity,  <i>  1*  0,  would  be  a  further  structure  above  and  beyond  these. 

It  is,  however ,  encouraging  to  see  that,  at  least  in  the  case  of 

linear  acceleration  the  formulae  (2.6)  and  (2.12)  can  accommodate  an 

inertial  field  ^  =  1 1>,  Z2i>  =  0  so  that  the  linear  acceleration 

can  be  conceived  as  a  gauge  transformation,  over  the  presence  of 

gravity,  as  it  was  always  assumed  it  would  be  the  case.  Incidentally , 

the  locality  condition  (5.3)  can  be  viewed  either  as  a  first  order 

statement  or  as  in  (2.6')  valid  in  second  order  as  well.  The  dif- 

ex  X  6  oi 

ference  contributes  to  L,  a  divergence,  %  J  =23  (<J>„3P4>, ) ,  hence  has 

r  a  a  p  A 

no  effect  on  the  theory. 

7.  QUANTIZATION 

The  classical  field  theory  of  gravitation  so  constructed  may  be 
quantized  by  simply  replacing  the  Lagrangian  (2.9)  with  an  operator 
expression 

Lf=:~  3^  3X$“  +  |  3X4>  3X*  .•  (7.1) 

where  semicolons  imply  normal  ordering.  By  a  well  known  procedure  of 
introducing  a  state  vector,  |  >,  a  unitary  transformation,  U,  etc., 
the  commutation  relations  can  lccally  be  found  in  the  usual  manner 
because  the  theory  is  locally  Lorentz  invariant.  In  other  words  the 
problem  of  quantization  of  a  gravitational  field  is  here  reduced  to 
the  quantization  of  a  special  relativistic  tensor  field,  i.  fact  an 
interesting  economy  in  formulation  is  realized  as  follows:  In  the 

usual  special  relativistic  case  two  independent  constraints  <  ]  -  4>  ]  > 

l  l  ^ 

and  <PV4)1J|!>  -  0  are  imposed  on  the  expectation  values.  First  of 

these  is  necessary  for  the  consistency  of  the  Lagrangian  procedure 

whereas  the  second  is  needed  for  the  consistency  of  quantization.  In 

the  new  theory  the  two  turn  out  to  be  equivalent  hence  only  one  is 
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necessary  to  impose.  The  reason  for  this  is  that  if  the  solution 

(2.6)  is  substituted  into  (2.1)  the  d'Alembert  operator  □ 2  is  the 

general  d ' Alembertian  whereas  when  we  contract  (2.1)  the  d'Alembertian 

5  y 

must  reduce  to  a  scalar  d'Alembertian  -is  to  be  interpreted 

as  a  scalar,  The  substitution  mentioned  yields 


2*V>  - 


trace  (  □  *$^) 


□J* 


4&  <|>d  <p 
v  a 


(7.2) 


v  y 

hence  3^$^  =  0  implies  ^  =  tj>  =  scalar.  In  the  process  of  quantization 

we  shall  therefore  need  only  the  condition ® 


<  ]  3  $V|> 

1  v  y  1 


=  0 


(7.3) 


The  commutation  relations  are 


$(x)  -  <t>(x)  l*(x')  =  ii(x-x')  (7.4) 

ik(x  -  x') 

D(x-x')  =  — —  f  dhke .  - -  (7.5) 

(2i r;4  c  k  -*2 

which,  by  a  suitable  substitution,  are  fo^nd  to  be  equivalent  to  a 
Gupta  type  of  quantization .9  Both  the  treatment  of  the  subsidiary 

conditions  (7.3)  and  the  deiivation  of  the  equations  of  motion  follow 
the  Gupta  procedure  with  similar  end  results.  The  process  of  quanti¬ 
zation  being  fairly  stbndard  the  main  improvement  in  the  subject  of 
quantized  gravity  is  here  due  to  the  improvement  of  the  background 
theory.  Since  there  exists  hardly  any  possibility  of  testing  the 
quantized  theory  in  the  foreseable  future  we  shall  emphasize  below 
only  some  formal  matters  which  seem  to  show  that  the  present  theory 
meshes  quite  naturally  with  the  microscopic  laws  of  quantum  mechanics. 

By  virtue  of  the  local  Lorentz  invariance  of  the  background  theory 
the  quantization  is  here  a  literal  adaptation  of  theusual  method  of 


quantization  with  the  added  feature  that  g^ (*,)  is  now  one  of  the  ob¬ 
servables.  The  metric  is  therefore  given  l y  its  expectation  value 


g  =  <|g  |>  =  <|n  •  e2<*  -  W\> 

yv  1 *yv 1  1  1 


(7.6) 


Conceptually  this  is  a  large  departure  from  the  Einoteinian  viewpoint 

where  g^  (x)  is  itself  the  gravitational  field.  The  following  comments 

might  be  of  interest:  a)  Defining  the  vacuum  by  $]>  -  0  and  assuming 

no  quanta,  real  or  virtual,  are  present  anywhere  the  expectation  value 

of  (7.6)  reduces  to  n^v  alone.  Hence  a  completely  empty  space  has  the 

property  of  global  Lorentz  structure,  b)  If  we  introduce  a  prescribed 

v 

classical  source  au^u  allow  a  large  number  of  quanta  (real  and/or 

virtual)  the  expectation  value  of  g ■  approaches  the  solutions  (time- 

dependent  and/or  static )  solutions  of  the  theory,  c)  Due  to  the  behavior 

-24> 

of  gQQ  as  a  negative  exponential  gQQ  -*•  e  the  theory  does  not  seem  to 
have  the  self  energy  difficulty  although  a  renormalization  might  still 
be  required  as  a  redefinition  of  the  rest-mass  or  of  g d)  In  prin¬ 
ciple  a  massive  gravitational  field  with  a  rest  mass  y  is  possible  to 
conceive  and  this  would  seem  to  obviate  the  necessity  of  dealing  with 
an  infrared  divergence  due  to  zero  graviton  rest-mass.  It  is,  however, 
difficult  to  justify  such  an  assumption  physically  as  gravity  seems  to 
be  such  a  long-range  force  that  it  applies  at  distances  of  intergalactic 

size,  e)  Since  in  the  absence  of  a  gravitational  stress-energy  tensor, 

v  v  v 

t  ,  the  two  formulae  D  (Ou  u  )  =  0  and  D  (au  )  -  0  are  equivalent  the 
y  v  y  v 

microscopic  theory  of  gravity  presented  in  this  section  may  be  conceived 
as  a  kind  of  spontaneous  symmetry  breaking,  f)  Quantized  form  agrees 
with  the  cross-section  for  gravitational  radiation  stated  previously^ . 
According  to  the  theory  gravitons  act  as  normal  spin-2  particles  with 
vanishing  rest  mass  and  standard  (i y^)  properties  of  polarization10 . 
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